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The  effects  of  charged  oxide  traps  on  the  electron  effective  mobility  (or 

conductivity  mobility)  and  the  linear  and  saturation  currents  in  the  inversion  surface 

channel  of  the  n-channel  metal-oxide-semiconductor  transistor  (nMOST)  are 

investigated.  The  deep  Ec-7eV  and  the  shallow  E^-leV  traps  of  the  oxygen  vacancy 

center  in  the  gate  oxide  are  charged  using  substrate-hot-electron  injection  (SHEi)  and 

Fowler-Nordheim  tunneling  electron  injection  (FNTEi)  methods  and  subsequently 

neutralized  by  electric  field-stimulated  electron  (EFSE)  emission  or  low-field  SHEi. 

The  effective  mobility,  fi^ff,  and  saturation  mobility,  n,^,,  of  the  electrons  are  monitored 

during  charging  and  discharging  cycles  and  shown  to  decrease  because  of  increased 

Coulombic  scattering  in  the  surface  channel  during  the  charging  stage,  and  to  recover 

during  the  discharging  stage.   In  addition,  changes  of  the  drain  saturation  current, 

^^D-sat'  ^i^d  gate  threshold  voltage,  AV^y,  are  also  monitored  and  the  effects  of  AVgt 


and  A|j.5^,  on  AI^.^^,  are  separated.  A  linear  dependence  of  the  electron  effective 
mobility,  Afi„^ff,  on  the  separated  positive  and  negative  charges  stored  in  the  oxide  is 
demonstrated  for  the  first  time. 

The  effects  of  areally  nonuniform  oxide  charges  and  interface  traps  in  the  gate 
oxide  are  separated  using  a  new  procedure  which  combines  two  experimental  d.c. 
characteristics:  the  subthreshold  drain-current  swing  and  the  d.c.  base  recombination 
current,  Ig.  Electrons  are  injected  into  the  gate  oxide  by  areally  uniform  and 
nonuniform  methods  to  demonstrate  this  new  technique. 

Finally  a  new  current-accelerated  channel-hot-carrier  (CHC)  methodology  is 
demonstrated  for  rapid  time-to-failure  extrapolation  of  submicron  MOSTs.  This  is 
demonstrated  using  the  n-channel  bipolar  MOST  (nBiMOST)  test  structure  with  Alg  as 
the  monitor  for  the  degradation.  This  technique  gives  an  acceleration,  or  stress  time 
reduction  factor  as  high  as  200  at  the  low  operation  drain  and  gate  bias  voltages  (less 
than  3.3V)  required  in  future  generations  of  deep  sub-half-micron  MOSTs.  This  is  in 
contrast  to  voltage-accelerated  channel-hot-electron  CHE  which  uses  higher  stress 
voltages  (6V,  8V  or  higher)  to  elevate  the  kinetic  energy  distribution  of  the  channel 
electrons  and  therefore  changes  the  dominant  failure  mechanisms. 


VI 


CHAPTER  1 
INTRODUCTION 


Since  the  initial  theoretical  work  on  the  carrier  mobility  in  the  semiconductor 
surface  inversion  layer  by  Schrieffer  [1,2]  and  the  channel  conductance  measurements 
of  Brown  and  Kingston  [3,4],  mobility  studies  have  been  made  extensively  using  the 
metal-oxide-semiconductor  field-effect  transistor  (MOST  or  MOSFET). 
Characterization  of  the  carrier  mobility  is  essential  for  advancing  today's  deep 
submicron  ULSI  (Ultra  large  scale  integrated)  circuit  technology  to  manufacture  multi- 
million-transistor  chips  since  it  is  directly  proportional  to  the  MOST's  output  current  to 
charge  up  the  load  capacitance  which  limits  the  clock  frequency  of  digital  ULSI 
circuits  and  since  it  is  also  directly  proportional  to  the  MOST's  transconductance  and 
the  cutoff  frequency  which  limits  the  bandwidth  of  analog  ULSI  circuits.  Accurate 
wide-range  characterization  of  the  carrier  mobility  requires  an  in-depth  quantitative 
understanding  of  the  physical  scattering  mechanisms  which  limit  the  mobility  in  the 
MOST  channel.  The  physics-based  understanding  is  essential  to  the  development  of 
accurate  MOST  device  simulators  that  do  not  rely  solely  on  calibrations  using 
experimental  data  that  cannot  assure  accurate  mobility  extrapolation  outside  of  the 
mobility  calibration  range. 

As  the  complementary-MOS  (CMOS)  and  bipolar-CMOS  (BiCMOS)  integrated 
circuit  technologies  continue  to  scale  down  the  gate  oxide  thickness  (below  50 A),  and 
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the  channel  lengths  (below  O.l^im)  [5-8],  the  impurity  dopant  concentrations  in  the 
CMOS  wells  and  the  MOST  channels  have  to  be  increased  accordingly  to  minimize  the 
various  short-channel  effects  which  give  undesirable  MOST  characteristics.  The 
electric  field  also  increases  as  the  dimension  decreases.  The  higher  electric  field  and 
dopant  impurity  concentration  confine  the  inversion  layer  charge  closer  to  the 
interface,  causing  a  higher  carrier  scattering  rate  due  to  the  proximity  of  the  carriers  to 
the  atomically  rough  oxide/silicon  interface  and  the  charges  in  the  gate  oxide.  The 
higher  electric  field  also  increases  the  electric  instability  of  the  gate  oxide  via  electron 
and  hole  trapping  at  the  existing  and  stress-generated  oxide  traps  and  at  the  stress- 
generated  oxide/silicon  interface  traps.  Furthermore,  the  higher  transistor  density  from 
smaller  transistor  dimensions  would  raise  the  operation  temperature  of  the  transistors 
and  the  ULSI  circuit  chips.  Thus,  the  three  major  mobility  reduction  factors  are  the 
higher  longitudinal  and  transverse  electric  fields,  the  confinement  of  the  charge  carriers 
(electrons  or  holes)  near  the  oxide/silicon  interface,  and  higher  operation  temperatures. 
The  effects  of  trapped  charges  in  the  silicon  substrate,  gate  oxide  layer,  and 
polysilicon  (poly)  gate,  and  at  the  poly-gate/Si02  and  Si02/Si-substrate  interfaces  on 
MOST  characteristics  have  been  studied  by  numerous  investigators  and  reported  in  the 
literature  [9].  The  magnitude,  sign,  and  distribution  of  these  trapped  charges,  in 
relation  to  the  conduction  carriers  (electrons  or  holes)  in  the  surface  channel,  can  alter 
the  threshold  voltage,  Vq^,  and  the  channel  effective  mobility,  both  of  which  govern 
the  shape  of  MOST's  drain  current-voltage  (Id-V^)  characteristics.  So  far,  studies  on 
the  effects  of  these  charges  on  the  carrier  mobility  have  been  limited  to  carrier 
scattering  by  the  charged  Si/Si02  interface  traps  and  ionized  impurity  centers  in  the  Si 


substrate,  and  the  Coulombic  screening  of  the  charged  traps  by  the  large  carrier  density 
in  the  surface  inversion  channel. 

In  general,  confinement  of  the  inversion  carriers  closer  to  the  interface  will 
increase  surface  (atomic)  roughness  scattering  and  convert  phonon  scattering  from  the 
three-dimensional  (3-d)  to  the  more  effective  2-d  'surfons'  scattering  as  the  quantum 
energy  states  in  the  potential  well  of  the  surface  channel  become  quantized  [10,  p.668]. 
Furthermore,  there  will  be  increased  Coulombic  scattering  from  the  charges  in  the 
poly-gate  and  the  oxide  and  interfacial  layers  mentioned  above  because  of  the  closer 
proximity  of  these  charges  to  the  carriers  in  the  conducting  channel.  Carrier 
confinement  will  reduce  scattering  from  ionized  bulk  impurities  unless  there  is  a 
corresponding  increase  in  the  impurity  concentration;  higher  impurity  concentration  is 
frequently  used  to  suppress  short  channel  effects  in  recent  generations  of  submicron 
MOSTs.  The  higher  MOST  operation  temperature  from  higher  dc  power  density  will 
also  increase  phonon  scattering  and  reduce  the  coulombic  scattering. 

There  have  been  two  effective  mobilities  measured  and  theorized  for  the  surface 
channel  and  the  MOS  transistors.  The  first  effective  mobility  is  obtained  from  the 
MOST's  channel  conductance  measurements  in  the  linear  channel  current- voltage 
range,  given  by  the  drain  conductance,  g^^dl^JdYj^  with  VQ=constant,  in  the  linear 
MOST  drain  current-voltage  (Id- Vq)  range.  This  will  be  denoted  by  pi^,,  li^near  ^r  |iiin- 

The  second  effective  mobility  is  obtained  from  the  MOST's  transconductance 
measurements,  g^=dlY)/dYQ  with  V^^constant,  which  will  be  denoted  by  [i^^^  This 
experimental  mobility  is  highly  dependent  on  the  energy  and  spatial  distribution  of  the 
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interface  traps  at  the  oxide/silicon  interface.  Thus,  it  is  difficult  to  obtain  reliable 
results  from  an  analysis  of  the  experimental  data. 

The  third  effective  mobility  is  obtained  from  least-squares-fit  of  the  experimental 
saturation  drain  current  versus  gate  voltage  data  (lo-sat'^c)  ^o  the  analytical  equation 
based  on  the  simple  bulk-charge  theory.  This  mobility  will  be  denoted  by  |J,^,. 

Two  effective  mobilities  are  measured  and  compared  with  theory  in  this  thesis: 
the  linear  and  saturation  mobilities,  |i,i„  and  [lg^^.  The  linear  mobility  provides  the  most 
direct  experimental  evaluation  of  the  various  scattering  mechanisms  which  reduce  the 
mobility.  The  saturation  mobility  gives  the  transistor  parameter  that  goes  into  the 
transistor  model  used  in  the  circuit  analysis  of  MOS  ULSI  circuits. 

In  this  thesis  the  effects  of  charge  trapping  and  detrapping,  in  the  gate  oxide 
layer,  on  the  two  effective  mobilities  are  studied.  Variation  of  the  drain  current  due  to 
controlled  changes  of  oxide  charge  density  will  be  correlated  individually  to  AY^-j-  and 
Afij^,.  This  gives  the  effects  of  Coulombic  scattering  of  the  channel  electrons  by  the 
positive  and  negative  charges  in  the  gate  oxide  on  mobility.  All  measurements  were 
undertaken  using  the  n-channel  bipolar-MOST  (n-BiMOST)  test  structure,  therefore, 
the  scope  of  this  thesis  is  limited  to  the  study  of  electron  effective  mobilities  in  the 
surface  channel.  Nevertheless,  the  new  fundamental  understanding  and  new 
experimental  methodologies  are  also  applicable  to  hole  mobilities  in  silicon  p-channel 
MOSTs. 

Results  of  detailed  studies  are  described  in  this  thesis  on  the  variation  of  the 
mobilities  {[i^^^  and  fij^j),  gate  threshold  voltage  V^-j-,  and  drain  saturation  current,  I^.^ap 
caused  by  the  three  charge  states  (-1,  0,  +1)  of  the  bridging  oxygen  vacancy  in  the  gate 


oxide.  The  oxygen  vacancy  was  also  known  as  the  E'  center  delineated  prior  to  1979 
by  electron  paramagnetic  resonance  experiments.  The  oxygen  vacancy  was  charged 
and  discharged  or  neutralized  by  Fowler-Nordheim  tunneling  electron  injection 
(FNTEi),  substrate  hot  electron  injection  (SHEi),  and  Fowler-Nordheim  tunnelling 
electron  emission  (FNTEe)  techniques,  while  the  MOST's  Iq-Vq  and  Id-V^ 
characteristics  are  monitored.  Koomen's  [11]  split-CV  method  was  used  to  measure 
the  carrier  capacitance  in  the  inversion  surface  channel,  Cj^y,  which  was  then  used  to 
calculate  the  inversion  carrier  charge  density,  Qjnv  The  linear  mobility  was  then 
computed  from  |i,i^  =  gy(W/L)Qj(^v  where  W/L  is  the  channel  width  to  length  ratio  or 
the  aspect  ratio.  The  linear  Iq-Vg  (at  three  constant  Y^:  V^j,  Vd2,  and  V^g)  and  the 
saturation  iD-^at-Vo-sa.  (Vq  =  Vg  =  Vp„^,,.,„pp,y  >  Vg  -  Vgt  =  Vp.,,,)  were  measured 
before  and  after  each  electron  injection-capture/emission  cycle  which  gave  the  Vgt, 
^D-sat'  M^sat  ^^'^  ^^Mn  ^^  ^  function  of  the  charged  oxygen  vacancy  center. 

The  second  part  of  this  research  concerned  the  effects  of  nonuniform  charging  of 
the  gate  oxide  using  channel-hot-electron  (CHE)  and  substrate-hot-electron  (SHE) 
injection  to  stress  the  n-channel  MOS  transistor.  A  new  technique  is  demonstrated  for 
separating  the  effects  of  uniform  or  nonuniform  interface  traps  and  nonuniform  oxide 
traps  along  the  length  of  the  MOST  channel.  The  technique  makes  use  of  the  recently 
developed  direct-current  current- voltage  (DCIV)  method  [12]  to  monitor  interface  trap 
generation,  and  the  well-known  change  of  the  subthreshold  drain-current  swing,  AS,  to 
monitor  the  spatial  uniformity  of  the  trapped  oxide  charge  during  the  controlled  oxide 
trap  (oxygen  vacancy)  charging  and  discharging  cycles.  In  addition,  the  recently 
demonstrated  current-acceleration  methodology  for  rapid  time-to-failure  measurement 


in  bipolar  junction  transistors  (BJT),  proposed  by  Neugroschel  and  Sah  [13],  is 
extended  for  the  first  time  to  MOS  transistors  under  CHE  stress.  This  new  MOST 
acceleration  stress  technique  makes  use  of  the  BiMOST's  unique  vertical  BJT  bottom- 
emitter  junction  to  accelerate  the  time-to-failure  while  maintaining  the  same  drain  and 
gate  voltages  (V^.  Vq).  which  determine  the  kinetic  energy  and  distribution  of  the 
injected  carriers.  This  new  MOST  acceleration  stress  technique  is  further  extended  to 
regular  MOSTs  which  do  not  have  a  substrate  junction  emitter  by  forward  biasing  the 
source  junction  as  proposed  by  Sah.  The  three  failure  parameters  measured  are  the 
increasing  base  current  measured  by  the  DCIV  method,  the  change  in  the  gate 
threshold  voltage,  and  the  change  in  the  drain  saturation  current. 

Chapter  2  will  review  the  scattering  mechanisms  affecting  the  carrier  mobility  in 
the  surface  inversion  layer  of  MOS  transistors.  The  measurement  setup  and  BiMOS 
transistor  fabrication  and  operation  details  will  be  described  in  chapter  3.  The  effects 
of  charged  oxide  traps  on  the  mobility,  saturation  current,  and  threshold  voltage  are 
discussed  in  chapter  4.  Experiments  demonstrating  a  new  measurement  technique  for 
the  separation  of  interface  traps  and  areal  oxide  charge  nonuniformity  will  be  presented 
in  chapter  5.  Finally,  the  two  new  current-accelerated  channel  hot  electron  stress 
methodologies  are  presented  in  chapter  6  which  forward  bias  the  substrate  or  source 
n/p  junction  emitter  to  increase  the  channel  hot  electron  current. 


CHAPTER  2 
SURFACE  INVERSION  CARRIER  TRANSPORT  THEORY 


The  three  primary  scattering  mechanisms  which  contribute  to  the  inversion 
carrier  mobility  degradation  as  reviewed  by  Sah  and  his  graduate  students  [14-16] 
included  the  following:  1)  lattice  scattering  due  to  acoustical  and  optical  phonons,  2) 
Coulombic  scattering  due  to  ions  and  dipoles  in  the  oxide  and  at  the  Si02/Si  interface, 
and  impurities  in  the  silicon  surface  space-charge  layer,  and  3)  surface  roughness 
scattering  due  to  the  interfacial  atomic  irregularities.  One  must  also  consider 
Coulombic  screening  by  the  inversion-layer  carriers  and  the  effective  mass  anisotropy 
in  2-d  conduction  due  to  1-d  quantization  in  the  surface  channel's  quantum  well  from 
the  high  transverse  electric  field.  These  scattering  mechanisms  will  be  reviewed  in  this 
chapter  based  on  the  theoretical  analyses  given  by  the  previous  investigators  as 
reviewed  by  Sah  [14-16]. 

2. 1  Theoretical  Calculation  of  li 

The  general  theory  of  inversion  layer  transport  and  carrier  effective  mobility  was 
first  presented  by  Schrieffer  [1]  to  account  for  the  reduced  conductivity  due  to  diffuse 
surface  scattering.  Diffuse  scattering  refers  to  the  assumption  that  the  probability  of 
scattering  is  equal  in  every  direction.   Schrieffer  solves  the  Boltzmann  transport 
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equation  using  the  relaxation  time  approximation  to  obtain  the  inversion  carrier 
distribution  function, 

f  =  fo  +  fi(v.z)  (2.1) 

where  fj  represents  a  small  perturbation  to  the  equilibrium  Boltzmann  distribution 
function,  fg.  Using  a  random  scattering  boundary  condition  to  eliminate  f,  at  the 
interface,  (x=0),  the  current  density  is  found  by  integrating  f  j  in  the  velocity  space, 

i^,  =  qjdv^dvydv,  Vyf,  (2.2) 

Here  v^  and  v^  are  the  carrier  velocities  parallel  to  the  interface  and  x  is  the  normal  (or 
transverse  direction)  to  the  interface.  By  integrating  the  current  density  from  the 
interface  to  the  edge  of  the  inversion  layer,  X;.  one  finds  the  total  d.c.  current,  ly,  due  to 
a  constant  longitudinal  field  Ey.  The  conductivity  effective  mobility  or  linear  mobility 
is  then  defined  as, 

h  =  qiiunNiKvEy.  (2.3) 

The  inversion  carrier  density,  Np^y,  is  calculated  for  an  assumed  triangular  potential 
well  (constant  transverse  field.  Ex), 

Ej,  -  dO/dx  =  constant  (2.4) 

V^O  =  p/e^  (2.5) 

or  by  solving  the  Poisson  equation  for  the  potential  well  profile  at  the  silicon  surface. 

2.2  Surface  Quantization 

Schrieffer  at  the  suggestion  of  Bardeen  [2]  also  recognized  the  possibility  of  1-d 
or  transverse  quantization  at  high  transverse  fields  E^  near  the  silicon  surface.  As  the 


applied  d.c.  gate  voltage  increases  the  surface  energy  band  bending,  the  potential  well 
at  the  silicon  surface  becomes  narrower  and  deeper;  eventually  quantizes  the  motion  in 
the  x-direction  giving  quantized  (discrete)  energy  subbands  perpendicular  to  the 
surface.  This  restricts  inversion  carrier  motion  to  the  two  dimensions,  y  and  z,  parallel 
to  the  interface  within  each  2-d  energy  subband.  The  two-dimensional  transport  in  the 
surface  inversion  layer  was  first  observed  by  Fang  and  Fowler  [17,18]  using 
magnetoconductance  measurements  on  a  <100>  silicon  surface  at  low  temperatures. 
Stern,  Howard  and  Fang  [17-20]  performed  selfconsistent  calculations  of  the 
Schroedinger  and  Poisson  equations  to  illustrate  the  energy  splitting  between  the 
ground  and  the  first  excited  2-d  energy  subbands.  Calculations  were  made  using  the 
variational  approach  with  the  assumption  that  all  the  carriers  occupy  the  lowest  energy 
subband  and  that  the  potential  is  a  function  of  x  only.  The  boundary  conditions 
assumed  were  as  follows:  (i)  the  potential  vanishes  at  x=oo,  (ii)  the  transverse  field, 
d^/dx,  vanishes  at  the  interface  (x=0),  and  (iii)  the  electron  wave  function  4^  vanishes 
at  the  interface  (infinite  SiOj/Si  barrier  height  or  zero  penetration  into  the  oxide)  and 
deep  in  the  bulk  silicon  (x=c>o).  The  envelope  function,  ^(z),  of  the  electron  wave 
function 

"¥  =  ^(z)exp(ikix  +  ik2y)  (2.6) 

is  then 

^(z)  =  C(z)exp[-iz(Wi3k,  +  W23k2)AV33]  (2.7) 

where  ^(z)  are  the  eigenfunction  solutions  of  the  Schroedinger  equation  and  W--  is  the 
reciprocal  effective  mass  tensor  (1-x,  2-y,  3-z).  Solving  the  wave  equation,  they 
compute  the  energy  eigenvalues  and  an  average  distance  from  the  interface  for 
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inversion  carriers  in  the  channel.  Their  results  illustrate  the  dependence  of  the 
quantized  energy  splitting  on  well  doping  concentration,  temperature,  transverse  field 
in  the  silicon,  E^,  and  crystal  orientation. 

Subsequent  theoretical  work  by  Pierret  and  Sah  [21]  supported  the  strong 
variation  in  effective  mobility  with  well  doping  concentration  in  weak  inversion  near 
Vg-p  (when  the  surface  potential,  O^,  equals  twice  the  bulk  potential  <I)5=2<J)p).  They 
also  noticed  that  the  effective  mobility  approached  a  common  value  versus  doping  at 
higher  transverse  fields  or  as  the  inversion  charge  density  became  degenerate. 

2.3  Crystalographic  Anisotropy 

In  the  surface  inversion  layer  the  silicon  crystal's  cubic  symmetry  is  disrupted  by 
the  transition  into  the  SiOj  layer  and  hence  one  must  also  refine  the  isotropic  bulk 
mobility  theory  to  include  the  crystal  orientation  dependence.  The  first  theoretical 
treatments  were  presented  by  Ham  and  Mattis  [22]  assuming  a  constant  surface  field 
and  using  a  special  set  of  spatial  transforms  to  convert  the  ellipsoidal  energy  surfaces 
to  spherical  ones.  Their  constant  transverse  field,  E^,  (or  triangular  potential  well) 
assumption  was  removed  by  Pierret  and  Sah  [21]  who  gave  the  solutions  for  the  exact 
electric  field  which  used  the  Boltzmann  transport  equation  and  considered  the 
distribution  function,  f,  for  each  of  the  six  ellipsoidal  conduction  valleys  in  E-k  space. 
The  analysis  was  then  converted  from  k-space  to  velocity  space  (v-space)  using  the 
spherical  relation  for  the  velocity, 

v^  =  (l/h)dE/dk  (2.8) 
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while  the  appropriate  mass  tensor  is  substituted  for  the  spherical  effective  mass,  m*. 
Pierret  and  Sah  also  assumed  that  the  surface  inversion  layer  is  nondegenerate,  there  is 
no  intervalley  scattering,  and  a  constant  relaxation  rime  may  be  used.  The  key 
difference  in  this  solution  compared  to  the  isotropic  case  is  that  the  inversion  carrier 
velocity  in  the  longitudinal  direction,  Vy,  is  now  dependent  on  the  transverse  velocity, 
v^.  Fortunately,  however,  if  the  normal  to  the  interface  is  in  the  <100>  direction,  v^,  is 
once  again  decoupled  from  the  carrier  velocity  parallel  to  the  interface  and  conduction 
becomes  isotropic.  Since  the  samples  tested  in  the  present  study  are  oriented  in  the 
<100>  direction,  as  is  the  case  in  most  of  today's  and  future  MOST  technologies,  the 
effects  from  anisotropic  conduction  will  not  be  included  in  this  thesis. 

2.4  Screened  Coulombic  Scattering 

As  discussed  earlier,  the  sources  of  Coulombic  scattering  include  1)  ionized 
impurities  in  the  surface  space  charge  layer  of  the  silicon  substrate  and  the  polysilicon 
gate,  2)  charges  in  the  gate  oxide,  and  3)  isolated  or  dipole  type  charges  located  at  the 
polysilicon-  gate/oxide  and  oxide/silicon-substrate  interfaces.  The  key  parameters 
governing  this  form  of  scattering  are  1)  Coulombic  screening  by  the  mobile  inversion 
charge  in  the  channel,  2)  the  distribution  of  these  mobile  charges  within  the  surface 
space-charge  layer  in  the  silicon,  3)  the  proximity  or  distribution  of  all  the  Coulombic 
scattering  sources  with  respect  to  the  inversion  layer,  and  4)  the  position  correlation  of 
the  scatterers  or  the  degree  of  randomness  in  their  spatial  distribution.  The  general 
problem  for  calculating  mobility  involves  finding  the  change  in  potential  energy  of  the 
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inversion  carriers  in  the  presence  of  these  scattering  centers.  For  example,  if  a  scatterer 
lowers  the  potential  energy,  then  the  carrier  density  increases  in  the  vicinity  of  this 
scatterer  while  the  reverse  is  true  for  a  decreased  potential.  Hence  the  mobile  charge 
redistributes  itself  so  that  the  scattering  center's  potential  is  screened. 

Greene  and  O'Donnel  [23]  first  considered  the  problem  of  scattering  in  the 
presence  of  surface  charges.  They  calculated  the  surface  reflectivity  in  terms  of  the 
differential  scattering  probability  for  a  shielded  Coulombic  potential  of  N  randomly 
placed  surface  charges.  Stem  and  Howard  [19]  also  reformulated  their  analysis  for  the 
lowest  energy  subband  in  the  quantized  surface  to  account  for  the  additional  screened 
Coulombic  potential  due  to  ionized  impurities  in  the  silicon  depletion  layer  as  well  as 
any  oxide  and  Si/Si02  interface  charges.  Using  first  order  perturbation  theory,  their 
assumed  potential,  *(r),  which  varies  with  r=(y^+z^)''''^  and  x,  induces  a  change  in  the  i* 
eigenvalue  5Ei(r)  of  the  mobile  inversion  charge, 

5Ei(r)  =  -qo(r)  =  -qJ5a)(r,x)g.(x)  dx  (2.9) 

g,(x)  =  (l/2)b\2gxp(-bx)  =  lC,(x)F  (2.10) 

b  =  {[487tq2m3/e,n][NDEp  +  (11/32)N,nv]}'^.  (2.11) 

where  gi(x)  and  ^^(x)  are  the  charge  distribution  function  and  the  normalized 
eigenfunction  for  the  Stern  and  Howard  2-d  ground  state  subband  [19].  The  induced 
charge  density  in  the  inversion  layer,  pj^j,  is  expressed  in  terms  of  the  reciprocal 
screening  length  Sj, 

Pi„,(r,x)  =  -  (ey27r)qli  s^  ^^(r)  gi(x).  (2.12) 

Ning  and  Sah  [24]  extend  this  single  ion  formulation  by  solving  Poisson's  equation  for 
SO  due  to  a  distribution  of  oxide  charge  at  r„  locations  in  a  plane  within  the  oxide. 
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V250(r,z)  -  2sO(r)g(z)  =  (47r/ejql5(r-r„)5(z)  (2. 1 3). 

NpxA  is  the  total  number  of  oxide  charges  in  a  sheet  of  area  A  with  an  oxide  charge 
areal  density  N„^.  Once  the  potential  fluctuation  0(r)  is  known,  the  scattering 
transition  rate,  T^^.,  is  found  using  the  'golden  rule'  and  the  Stem  and  Howard  [19]  2-d 
ground  state  plane  wave,  4^, 

r^.  =  27t/n  <l<4',,|(-q5<D)l4^,>|2>^^^.  (2.14) 

The  inversion  layer  mobility,  |I,  may  then  be  found  numerically  or  by  using  the  Bom 
approximation  to  calculate  F^., 

[i  =  (q/mjjyt(f„  -  l)dy/lf„dy  (2.15) 

T  =  Zrkk.(l-cos0).  (2.16) 

Stem  and  Howard  [19]  illustrate  the  use  of  the  Bom  approximation  for  a  single  ionized 
impurity  located  at  x=x^  and  discuss  the  validity  of  this  approximation  in  great  detail 
by  comparison  with  numerical  phase-shift  calculations.  The  approximation  is 
inaccurate  at  extremely  low  temperatures,  or  high  doping  impurity  concentrations  or 
other  scattering  charge  densities.  This  stems  from  the  first  order  perturbation  analysis 
used  in  this  approximation  which  assumes  the  2-d  wave  function,  4^,  is  unchanged  in 
the  presence  of  the  charged  scatterers.  Other  sources  of  error  found  in  this  analysis 
include  the  assumptions  of  a  spherical  effective  mass,  the  neglect  of  screening  by 
inversion  carriers  and  the  possibiUty  of  conduction  in  higher  energy  subbands  at  higher 
transverse  fields,  temperatures  and  inversion  carrier  densities.  Furthermore,  the 
average  distance  between  the  inversion  layer  and  the  surface  charge  scattering  centers 
was  shown  to  have  a  pronounced  effect  on  the  mobility  as  would  be  expected  from  the 
1/r  dependent  Coulombic  potential. 
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The  subject  of  the  charged  scatterer's  spatial  distribution  and  the  effects  of 
position  correlation  of  the  scatterer  were  treated  by  Ning  and  Sah  [24,  25].  These 
authors  reformulate  the  scattering  problem  into  a  simpler  2-d  case  where  the  potential 
perturbation  arises  from  the  spatial  variations  in  the  local  charge  density.  The  charge 
density  fluctuations  5p  are  represented  as, 

5p(r)  =  p(r)-p"(x)  (2.17) 

p°(x)  =  <p(r)>,,,  (2.18) 

where  p°(x)  represents  the  x  spatial  dependence  of  p(r)  averaged  over  the  y-z  plane.  If 
the  thermal  wavelength  of  the  inversion  carriers,  \y^,  is  much  larger  than  the  inversion 
layer  thickness,  X;,  then  the  inversion  layer  can  be  treated  as  a  2-d  gas.  Under  these 
conditions  the  perturbation  analysis  proceeds  as  before  to  solve  for  the  scattering 
transition  rate  and  mobility.  Three  sample  oxide  and  inversion  charge  distributions, 
with  or  without  the  effects  of  screening  and  correlation  were  analyzed  [24]  and  used  to 
analyze  the  data  [25,  26].  Their  results  indicate  that  increased  correlation  and/or 
screening  by  the  inversion  charge  would  increase  the  channel  mobility. 

The  carrier  scattering  by  neutral  dipoles  formed  by  the  hydrogen  passivation  of 
ionized  surface  states  at  the  Si/SiOj  interface  as  first  suggested  by  Sah  also  falls  under 
the  analysis  discussed  here.  This  dipole  scattering  rate  was  derived  by  Hess  and  Sah 
[27]  using  the  aforementioned  Born  approximation.  Their  results  indicate  a  T'^ 
dependence  in  the  mobility  for  T>100°K, 

Hdp  =  HJ'-V(NdpLdp2)  (2.19) 

where  Npp  is  the  density  of  dipoles  at  the  interface  and  Lpp  is  the  effective  spatial 
separation  of  the  positive  and  negative  charges  in  each  dipole. 
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2.5  Phonon  Scattering 

As  the  temperature  of  the  semiconductor  increases,  scattering  due  to  quantized 
lattice  vibrations,  or  phonons,  will  naturally  increase.  These  lattice  vibrations  are 
commonly  categorized  as  acoustical,  for  those  with  lower  energies,  and  optical  (or 
intervalley)  for  higher  energies.  The  first  theoretical  analysis  for  surface  acoustical 
phonon  scattering  by  Kawaji  [28]  was  essentially  a  2-d  version  of  the  Bardeen- 
Shockley  deformation  potential  theory  for  bulk  semiconductors.  In  this  analysis  a 
change  in  potential  (assumed  to  be  an  inverted  triangular  well)  at  the  Si/Si02  interface 
originates  from  longitudinal  phonons  whose  wave  vectors  are  confined  to  the  2-d  plane 
parallel  to  the  interface.  Hence  these  phonons  only  interact  with  inversion  carriers 
moving  parallel  to  the  interface.  The  lattice  displacement,  5R(R„),  in  a  direction  1^  and 
the  corresponding  dilation,  A(r),  due  to  a  lattice  wave  in  the  2-d  surface  plane  are  given 

by, 

5R(R„)  =  N^-"^  1,  [a,exp(ikR„)  +  a;exp(-ikR,)]  (2.20) 

A(r)  =  V-5R(R„).  (2.21) 

There  are  N2  vibrating  lattice  points,  and  L  represents  the  number  of  atomic  layers  in 
the  conducting  surface.  The  semiconductor  density,  p,  and  the  mass  per  unit  area,  M, 
in  the  surface  inversion  layer  of  thickness  'd'  are  then  related  as  p=MLN2.  The  2-d 
deformation  potential  is  given  by,  A(r)E2,  where  E2  is  the  deformation  potential 
constant  which  was  obtained  by  fitting  the  theory  to  experimental  mobility  data.  This 
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potential  may  now  be  used  to  calculate  the  scattering  transition  matrix  element,  Mpp-, 
between  the  quantized  2-d  inversion  carrier  momentum  states  p  and  p', 

Mpp.  =  J4^|,(p)A(r)E24^,|(p')da  (2.22) 

p  =  p'±k     (2.23) 

x^  =  !m\  1  -  cos0)V2d?i  (2.24) 

The  relaxation  time.  T,  refers  to  the  electrons  traveling  parallel  to  the  interface  with 
wave  functions  4^,1  and  4^ * ,  where  Vj  is  the  2-d  density  of  states  along  a  constant  energy 
curve,  A,,  in  E-k  space.  Assuming  a  constant  velocity,  Cjj.  for  the  2-d  phonon  wave  and 
an  effective  longitudinal  mass,  m,,,  for  the  inversion  carriers  in  the  2-d  energy  subband, 
the  mobility  is  calculated  as, 

^L  =  {q^WnmK  E^keT).  (2.25) 

Considering  only  the  2-d  ground  state  eigenvalue  E^(0)  for  the  triangular  potential 
well,  and  given  the  transverse  field  as. 

Ex  =  Ex(0)/qd  =  (47r/eJ(QDEP  +  Qinv)  (2.26) 

it  is  easily  shown  that  d  (the  average  distance  of  the  inversion  layer  charge  from  the 
interface)  has  a  (N^gp  -i-  Nj;^)""^  dependence  and  hence  the  mobility  is  proportional  to 

T^HNdep  +Ninv)-^^- 

The  theory  was  later  refined  by  Ezawa,  Kawaji  and  Nakamura  [29]  who 
introduced  the  term  'surfon'  to  describe  acoustic  phonons  that  satisfy  the  interface 
boundary  conditions,  but  have  bulk  and  surface  modes  of  vibration.  Their  analysis 
included  various  vibrations  modes  for  the  surtons  and  their  electron  interactions,  as 
well  as  the  effects  of  anisotropic  conductivity  tensors  for  different  crystal  orientations. 
Inter-subband  transitions  and  intervalley  scattering  were  not  considered.  Some  of  the 
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remaining  discrepancy  in  the  theory,  on  the  order  of  6X  larger  mobility  than  that 
measured,  was  attributed  to  optical  phonon  interactions  and  possible  differences 
between  the  assumed  bulk  values  of  the  constants,  such  as  polarizability,  e,  and  their 
actual  values  near  the  Si/Si02  interface. 

For  temperatures  below  100°K  the  measurements  of  Sah,  Ning  and  Tschopp  [25] 
showed  that  the  mobility  due  to  lattice  scattering,  ji^,  varies  as  7.4xlO^/T  by  using  the 
approximation  by  Debye  and  Conwell  for  |ij  [30], 

1/^lL  =  1/llexp  -  Wl  (2.27) 

where  [i  and  fij  are  the  experimental  and  ionized  impurity  values  of  mobility.  At 
higher  temperatures  (>150°K)  Sah-Ning-Tschopp  [25]  reformulated  the  theory  to 
account  for  optical  or  intervalley  phonon  scattering  and  found 

Hl=10VtI  (2.28) 

This  solution  was  obtained  by  replacing  the  deformation  potential,  the  phonon  energy, 
and  the  occupation  number  for  acoustical  phonons  with  their  optical  counterparts 
which  was  also  employed  by  Luong  and  Shaw  [31].  The  two  relaxation  times  for  the 
optical  and  acoustical  cases  are  assumed  independent, 

x-'=xl'+Xo\  (2.29) 

The  resulting  expression  for  mobility  is  given  in  terms  of  the  purely  acoustical 
dependence  11^=7. 4xlO^/T,  the  optical  and  acoustical  ratio  of  the  deformation  potential 
constants  Z^/Z^,  and  the  optical  phonon  energy,  HcOq.  A  good  fit  to  their  mobility  data 
was  obtained  for  T=20-300°K  and  a  range  of  oxide  and  interface  trap  densities  when, 

Zo/Za  =  2.3  (2.30) 

hcOq/kg  =  650±50°  K  (2.31) 
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2.6  Surface  Roughness  Scattering 

The  possibility  of  mobility  reduction  due  to  the  atomic  irregularity  at  the  Si/Si02 
interface  was  analyzed  by  Cheng  and  Sullivan  [32-34].  Under  very  controlled 
processing  conditions,  they  were  able  to  fabricate  transistors  using  a  variety  of 
oxidation  conditions  to  alter  the  quality  of  the  Si/Si02  interface  and  a  range  of  interface 
trapped  charge  densities,  Njj.  Measurements  were  performed  at  liquid  helium 
temperature  (4.2°K)  to  minimize  the  interference  of  phonon  scattering.  Mobility  was 
extracted  at  high  transverse  fields  or  high  inversion  carrier  densities  to  emphasize  the 
effects  of  surface  roughness  by  bringing  the  mobile  carriers  closer  to  the  interface.  The 
observed  strong  Nj,^  (-1  <  a  <  -1.5)  dependence  of  the  surface  mobility  at  high  fields 
was  attributed  by  these  authors  to  atomic  surface  inhomogeneity. 

The  analysis  for  modeling  this  scattering  mechanism  is  similar  to  that  of  surface 
charges  with  a  random  spatial  distribution  in  the  oxide.  In  an  ideal  situation  the 
Si/Si02  interface  would  be  treated  as  a  perfect  plane.  Differences  in  the  crystal 
structure  of  the  Si  substrate  and  the  amorphous  nature  SiOj,  are  immediately  evident  in 
their  energy  bands  as  well  as  their  electrical  and  thermal  conductivity  properties  which 
prevent  the  Si/Si02  interface  from  being  completely  planar.  To  model  the  apparent 
atomic  irregularity  in  the  direction  perpendicular  to  the  interface,  a  random  spatial 
displacement,  Ax(R),  is  imposed  upon  the  potential  V(x)  of  an  ideal  planar  surface. 
The  perturbed  Hamiltonian  is  given  by, 

H  =  -  h^V^Ilm  +  V[x  -  Ax(R)]  (2.32) 

5V  =  V[x  -  Ax(R)]  -  V(x)  =  Ax.0V/3x)  (2.33) 
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and  5V  is  the  first  order  perturbation  of  the  potential.  Assuming  the  same  form  for  the 
ground  state  wave  function,  H',  as  Stern  and  Howard  [19],  though  translating  the 
boundary  condition  of  y(x=0)=0,  the  matrix  element  for  5V,  <C(x)l5VI^(x)>  and  the 
scattering  rate,  T^^.,  may  be  calculated.  The  derived  relaxation  time,  t(k),  exhibits  a 
strong  dependence  on  the  inversion  carrier  concentration  and  two  adjustable 
parameters  related  to  the  Gaussian  distribution  assumed  for  Ax(R), 

l/x=   (1447t2q^L2A2m|/e/h^)[NDEP  +  d  1/32)Ninv]' 

*j(l  -  cose)exp[-  (1  -  cosQ)kW2]  (2.34) 

where  A  is  the  mean  asperity  of  the  perturbation  to  the  ideal  plane  and  L  is  the 
correlation  length  of  the  Gaussian  distribution  (Figure  6  in  [32]).  Assuming 
completely  independent  scattering  mechanisms,  the  reciprocal  relaxation  times  for 
surface  roughness  and  the  Coulombic  scattering  are  added  and  used  to  find  the 
mobility  in  the  usual  manner, 

l/x=l/t,,+  l/x,„„,  (2.35) 

^i  =  (q/m)<x>  (2.36) 

Since  various  surface  preparations  have  been  used  in  the  industry  and  reported  in  the 
literature,  and  since  the  thermal  oxide  growth  conditions  (dry,  wet,  with  or  without 
HCl,  TCA,  and  nitrogen,  varying  temperature,  etc.)  also  vary  from  one  process  to  the 
next,  large  variations  in  the  asperity  and  correlation  would  be  expected.  Hence  it 
would  not  be  unusual  to  find  significant  process  variation  in  [i^  and  its  dependence  on 
NjNv.  Matsumoto  and  Uemura  [35]  refined  the  work  of  Cheng  et.  al.  [33]  by  including 
the  effects  of  screening  by  the  inversion  layer  and  by  selfconsistent  calculation  of  V(x) 
using  the  variational  approach  and  the  Stern  and  Howard  [19]  ground  state  wave 
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function.  They  also  defined  an  expression  for  the  effective  field  in  the  semiconductor 
as, 

E,^  =  JC(x)5VC(x)dx  =  4;rqVe,[(NiNv/2)+NDEp]  (2.37) 

where  the  first  term  (inversion  charge)  in  the  parentheses  will  dominate  at  high  fields 
as  expected.  The  relaxation  time  can  therefore  be  expressed  as  a  function  of  E^^. 
Harstein,  Ning,  and  Fowler  [36]  have  also  measured  the  effective  mobility  at  very  low 
temperatures  while  varying  the  oxide  charge  in  their  samples  and  have  experimentally 
separated  the  effects  of  Coulombic  and  surface  roughness  scattering  using  the 
Mathiesiens  rule.  More  recent  measurements  by  Takagi  [37]  at  77°K  have  shown  close 
agreement  to  the  theoretically  predicted  (Nj^v)"^  dependence  for  electron  mobility 
due  to  surface  roughness  whereas  the  mobility  for  holes  demonstrates  a  weaker 
(NiNv)"'  dependence.  Agostinelli  [38]  speculates  that  this  may  be  due  to  the  difference 
in  the  distribution  of  holes  in  the  inversion  layer  further  away  from  the  interface  as 
suggested  by  Moglestue's  [39]  self-consistent  calculations  on  <100>  siUcon. 

2.7  Physical  and  Semi-empirical  Models 

It  would  appear  that  out  of  the  three  different  scattering  mechanisms  affecting 
the  inversion  layer  mobility,  surface  roughness  scattering  is  the  most  difficult  to 
accurately  model  without  experimental  parameter  fitting  because  of  surface  processing 
dependences.  Yet  it  should  be  noted  that  some  of  the  recent  experimental  work  has 
shown  a  universal  behavior  for  surface  roughness  scattering  in  devices  processed  by 
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different  groups  [37,  40-42].  In  providing  a  physically  based  mobility  model,  Sah  [10, 
14-16]  summarized  the  various  scattering  mechanisms  affecting  mobility  as  follows: 
Electrons 


^iOnLA  =  7.4xl05n^ 
^OnLO=l-0xl0W-9 

^ionOX=10'(3xlO'VNox)(T/80) 

^iOnDP=10'(NDpLD2)-l(T/300)'-5 

^o„sr=1-5x102V/(IQbI+IQinvI)' 
^lonci  =  8.5x10V(Nion)''' 


Surface  acoustical  phonons  (2.38) 

Surface  optical  phonons  (2.39) 

Oxide  charge  (2.40) 

Surface  dipoles  (2.41) 

Surface  roughness  (2.42) 

Ionized  impurities  (2.43) 


where  Qg  is  the  bulk  charge  density  and  Nqx  is  the  oxide  charge  density.  The  total 
mobility,  [i^,  is  found  by  Mathiessen's  rule  by  adding  the  mobility  terms  listed  above 
as  |i''.  The  hole  mobilities  were  approximated  by  scaling  the  electron  mobility  by  the 
hole/electron  bulk  mobility  ratio  (Hp/fi^=478/1423)  [14-16].  The  transverse  and 
longitudinal  field  dependences  are  also  given  [10]  as  follows: 

EcL=0sat/^io  =  l-O2xlOX 

EcL=0sa/^io  =  O.75xlO>„ 

Ect=  100  kV/cm 
where  Ey  and  E^  are  the  longitudinal  and  peak  transverse  fields  in  the  silicon  surface 
layer,  E^l  and  E^^  are  the  critical  longitudinal  and  transverse  fields  and  the  constants,  y 
and  5,  are  empirical  constants  which  vary  slightly  with  temperature  and  both 
longitudinal  and  transverse  electric  fields. 


Longitudinal  Fields 

(2.43) 

Transverse  Fields 

(2.44) 

Electrons 

Holes 
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Various  other  semi-empirical  mobility  models  have  been  developed  for  electrons 
[43-49]  and  holes  [38]  but  as  the  term  semi-empirical  suggests,  these  models  require  a 
priori  knowledge  of  the  device  characteristics  to  establish  the  various  fitting  paiameters 
in  the  model.  The  model  developed  by  Schwarz  and  Russek  [43]  which  showed 
excellent  agreement  with  the  experimental  data  of  various  groups  [25,  40,  42,  50]  was 
extended  by  Shin  et  al.  [45,  46]  using  a  novel  modeling  approach,  to  account  for  the 
low-field  and  high-field  mobility  drop-off.  The  model  by  Shin  et  al.  model  extracts  the 
functional  dependence  of  the  inversion  layer  mobility  on  local  transverse  and 
longitudinal  electric  fields,  substrate  doping  concentration,  fixed  interface  charge  and 
temperature  from  the  experimentally  measured  effective  mobility,  fi^ff.  Their  model, 
which  was  implemented  in  the  PISCES  2-d  device  simulator,  was  shown  to  have  good 
agreement  with  experimental  Id-V^  and  Iq-Vq  data.  These  models  though  do  not 
account  for  the  Coulombic  scattering  due  to  trapped  oxide  charge  build-up  during 
stress  or  the  distribution  of  this  charge  which  are  the  primary  topics  of  this  thesis. 

The  following  chapter  will  the  discuss  the  theory  of  operation  and  fabrication  of 
the  BiMOST  test  structure  and  various  techniques  used  to  extract  the  saturation  and  the 
conductivity  effective  mobilities. 


CHAPTER  3 

BiMOST  FABRICATION  AND  OPERATION  AND 

EFFECTIVE  MOBILITY  EXTRACTION 


In  chapter  2  the  work  of  previous  authors  on  the  scattering  mechanisms  affecting 
the  inversion  carrier  effective  mobility  was  reviewed.  This  chapter  consists  of  two 
sections  on:  1)  the  fabrication  and  operation  of  the  BiMOST  test  structure  and  2)  a 
detailed  description  of  the  mobility  extraction  procedures.  The  following  chapter, 
chapter  4,  will  describe  their  applications  to  the  experiments  on  the  effects  of  charging 
and  discharging  oxide  traps  on  the  degradation  of  the  mobility,  drain  saturation  current, 
and  threshold  voltage. 

3.1  BiMOST  Fabrication  and  Operation 

The  BiMOST  test  structure,  shown  in  Fig.  3.1,  was  employed  in  this  study 
because  of  the  versatility  it  provides  for  studying  the  gate  oxide  layer.  The  physical 
aspects  which  make  this  structure  unique  are  that  it  combines  a  vertical  n/p/n  or  p/n/p 
bipolar  junction  transistor  (BIT)  with  an  n-  or  p-channel  MOST  on  the  surface.  This 
allows  one  to  stress  the  SiOj/Si  interfacial  layer,  inject  charges  into  the  gate  oxide  and 
characterize  the  resultant  degradation  by  methods  not  available  using  the  MOS  or 
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Figure  3.1    N-channel  BiMOST  test  structure  biased  for  substrate  hot 
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bipolar  transistors  or  the  MOS  capacitor  alone.  These  features  will  be  discussed  in  the 
following  section. 

The  BiMOST  structure  is  in  fact  present  in  all  CMOS  and  BiCMOS  technology 
today.  A  typical  CMOS  wafer  and  chip  fabrication  process  begins  with  a  highly  doped 
p-type  silicon  wafer  of  6-inch  or  8-inch  diameter  and  500-micron  thick  (to  be  denoted 
by  p+Si  where  +  signifies  high  dopant  impurity  concentration  which  is  the  boron 
acceptor  for  p-type  Si).  On  the  atomically  flat  and  polished  wafer  surface  is  epitaxially 
grown  a  thin  lowly  doped  p-type  silicon  layer  (to  be  denoted  by  p-Si  layer  where  the  - 
or  minus  sign  signifies  low  dopant  impurity  concentration).  Using  thermally  grown  or 
deposited  silicon  dioxide  (SiOj),  silicon  nitride  (Si3N4),  and  polycrystalline  silicon 
films  as  a  penetration  mask  and  photolithography  for  pattern  definition,  the  epitaxial 
surface  is  selectively  implanted  with  n-type  (phosphorus  or  arsenic)  and  p-type  (boron) 
impurities  to  form  the  n-well  and  p-well  in  which  the  p-channel  and  n-channel 
MOST's,  respectively,  are  then  built  respectively.  To  achieve  a  good  interfacial  layer 
for  oxide  growth,  improve  the  surface  mobiUty  and  minimize  short  channel  effects,  the 
impurity  profile  in  the  wells  is  highly  retrograded  in  today's  production  MOSTs  with 
<0.5micron  channel  length  and  in  future  state-of-the-art  0.1|im  channel  length  MOSTs 
[7].  A  typical  retrograde  profile  will  have  a  peak  impurity  concentration  on  the  order 
of  -lO'^cm'^  at  about  500A  from  the  Si/Si02  interface  in  the  ion  implanted  silicon  well 
while  the  impurity  concentration  at  the  Si/SiOj  interface  is  1-2  orders  of  magnitude 
below  this  peak.  If  the  inversion  layer  of  the  MOST  surface  channel  is  about  30-lOOA 
thick,  as  discussed  in  chapter  2,  then  the  highly  retrograde  profile  of  ionized  impurity 
centers  in  the  well  will  play  an  important  role  in  determining  the  effective  mobility. 
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Summarized  briefly,  the  following  steps  in  the  CMOS  (fabrication)  process  would  first 
thermally  grow  the  very  thin  (3.5nm  for  next  generation  and  6nm  for  Pentium  P-6 
announced  by  Intel  in  December  1994)  gate  oxide  on  the  epitaxial  surface  and  then 
deposit  a  polysilicon  layer  (poly-Si  or  just  poly)  over  the  entire  surface.  Next,  using 
the  poly  as  the  ion  implantation  mask  by  etching  holes  or  patterns  through  the  poly  and 
the  gate  oxide  layers  using  photolithography  (known  as  the  self-aligned  process) 
shallow  n-type  and  p-type  impurity  ion  implantation  through  the  holes  is  made  to  give 
the  n-i-  and  p-i-  source  and  drain  islands  and  to  dope  the  polysilicon  gate  to  high  n+  or 
p-i-  conductivity.  On  a  p-i-Si  or  n-i-Si  6"  or  8"  wafer,  both  a  p-channel  BiMOST  with  a 
vertical  p/n/p  BJT  and  an  n-channel  BiMOST  with  a  vertical  n/p/n  BJT  can  be 
fabricated.  A  cross-sectional  view  of  the  nBiMOST  is  shown  in  Fig. 3.1  which  is 
fabricated  on  a  n-i-Si  wafer. 

Some  of  the  MOSTs  and  BiMOSTs  measured  in  this  thesis  investigation  were 
obtained  from  industrial  mentors.  They  were  completely  processed  in  their  state  of  the 
art  production  facilities,  hence,  represent  the  cleanest  condition  attainable.  Large  area 
pBiMOSTs  and  nBiMOSTs  were  also  fabricated  by  this  candidate  in  our  university 
clean  room  fabrication  laboratory,  at  high  yields,  and  by  other  graduate  students  and 
professors  in  our  research  group.  These  large  area  transistors  were  fabricated  on  l"xl" 
wafers  which  were  cut  from  6"  or  8"  diameter  wafers  provided  by  industrial  mentors 
who  first  oxidized  the  n/p-epitaxy/p-i-substrate  and  p/n-epitaxy/n-i-substrate  6"-8" 
wafers.  All  the  remaining  fabrication  steps,  except  ion  implant,  were  carried  out  by 
this  candidate  in  our  own  clean-room  device  fabrication  laboratory,  which  included 
self-aligned  process  to  define  the  gate,  implant  the  source,  drain,  and  poly-gate  regions 
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(by  a  California  silicon  chip  foundry),  and  evaporate  aluminum  contacts  on  the  front 
and  back  of  the  wafers.  The  industrial  mentor  supplied  MOSTs  were  rectangular  with 
L=W=100|Lim  to  L  and  W  less  than  0.5|im.  Those  fabricated  by  this  candidate  were 
large  concentric  circles  with  L=150|im  and  a  large  gate  area  (3.89X10'^cm^).  Gate 
oxide  thickness  covered  the  range  of  3.5nm  to  30nm. 

One  of  the  key  advantages  of  the  BiMOST  test  structure  is  that  the  oxide  electric 
field,  Eqx'  the  reverse  bias  voltage,  V^g,  across  the  inversion  layer  or  the  collector/base 
junction  and  the  injection  current  through  the  gate,  1^,  can  all  be  varied  independently. 
This  allows  one  to  stress  the  oxide  using  the  substrate  hot  electron  injection  (SHEi) 
technique  at  gate  bias  voltages  as  low  as  V(3=Vqj  near  the  threshold  and  well  below  the 
Fowler  Nordheim  tunneling  voltage  or  tunneling  electric  field  .  For  SHEi  the  vertical 
n/p/n  bipolar  transistor  is  operated  in  the  forward  active  mode  of  the  bottom  emitter 
configuration  with  the  bottom  emitter-base  junction  (Vgg)  forward  biased  and  the  n- 
channel/p-well  space  charge  layer  reverse  biased.  The  forward  emitter-base  bias  (Vgg) 
injects  the  minority  carriers  (electrons)  from  n-emitter  into  the  p-base  of  the  vertical 
n/p/n  BJT.  These  minority  carriers  (electrons)  diffuse  through  the  p-base  and,  upon 
entering  the  surface  space  charge  layer  beneath  the  n-inversion  channel,  are  accelerated 
towards  the  Si/SiOj  interface  by  the  electric  field  in  the  surface  space-charge  layer  of 
the  reversed  biased  collector/base  and  channel/base  junction.  If  the  accelerated 
electrons  have  enough  kinetic  energy  at  the  interface  to  surmount  the  Si/Si02  electron 
energy  barrier  (3.13eV  for  electrons  and  4.25eV  for  holes),  they  will  be  injected  into 
the  oxide  conduction  band  (or  valence  band  for  holes  injected  in  the  pBiMOST).  In  the 
oxide  layer  the  injected  electrons,  now  labeled  oxide  electrons  or  electrons  in  the  oxide 
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conduction  band,  are  accelerated  further  by  the  oxide  electric  field,  Eq^,  as  they  transit 
toward  the  poly-Si  gate  electrode  and  exit  the  gate  terminal  to  give  the  measured  gate 
current.  Only  a  small  fraction  of  the  injected  oxide  electrons  will  be  trapped  or 
expended  via  one  of  numerous  trap  charging-discharging,  generation,  neutralization, 
and  even  annihilation  mechanisms  in  the  oxide  or  at  the  poly-gate/Si02  and  Si02/Si 
interfaces.  During  the  stress  or  SHEi,  the  gate  current  is  maintained  at  a  constant  value 
by  varying  the  d.c.  forward  bias  on  the  emitter/base  junction.  In  former  reports, 
verified  in  this  thesis,  SHEi  is  areally  uniform,  across  the  length  of  the  channel,  in 
contrast  to  the  ac  avalanche  injection  or  the  dc  channel  hot  carrier  (CHC)  or  channel 
hot  electron  (CHE)  stress  techniques.  Areal  uniformity  greatly  simplifies  the  analysis 
of  the  MOST  IV  and  CV  characteristics  since  any  distortion  of  the  IV  and  CV  curves 
due  to  nonuniform  charge  in  the  oxide  can  be  misinterpreted  as  interface  traps 
generated  during  the  stress.  By  maintaining  the  reverse  bias  voltage  on  the 
collector/base  junction  below  breakdown,  SHE  injection  remains  uniform  and  enables 
us  to  evaluate  the  initial  quality  (or  oxide  and  interface  trap  densities)  of  the  oxide  and 
Si/Si02  interface,  and  also  the  effects  of  post  oxidation  processing  undertaken  by  us. 

3.2  Measurement  of  u,i„ 

The  effective  mobilities  (linear  and  transconductance)  of  electrons  and  holes  has 
been  measured  extensively  in  the  literature  using  the  MOST.  Leistiko,  Grove  and  Sah 
[51]  measured  the  dc  or  differential  conductance,  g^,  of  the  inversion  layer  at  low  drain 
voltages  or  in  the  linear  I^-Vp  range,  which  was  defined  as 
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gd  =  (W/L)qja(x)dx  (3.1) 

where  g(x)  is  the  conductivity  of  the  inversion  layer  at  a  depth  x  below  the  interface 
defined  by  a(x)=q|i(x)N(x)  where  |i(x)  and  N(x)  are  the  electron  mobility  and  electron 
concentration  at  a  distance  x  from  the  Si02/Si  interface  (x=0).  Thus,  a(x)  is  the 
electron  conductivity  (siemen  or  mho/cm)  at  a  depth  x  in  the  inversion  channel.  The 
length  and  width  of  the  channel  are  L  and  W.  The  integration  covers  the  inversion  layer 
thickness  (x=0  to  Xj).  The  drain  conductance,  effective  conductivity  or  linear  mobility, 
and  inversion  charge  density  (q/cm^)  are  then, 

gd  =  (W/L)qj^(x)N(x)dx  =  (W/L)^l„„Qjf^  (3.2) 

^^Hn  =  gd/[(W/L)QiNv]  =  j^(x)N(x)dx  /  NjNv  (3.3) 

NiNV  =  iN(x)dx  =  QiNv/q.  (3.4) 

Approximating  L  and  W  of  the  large  area  devices  by  their  drawn  length  and  width,  one 
needs  only  to  measure  the  drain  conductance,  g^,  and  inversion  charge,  Qj^v,  to  extract 

In  this  investigation,  the  drain  conductance  is  calculated  from  the  experimental 
dc  Id-Vjjs  characteristics  in  the  linear  region  (V^s  =  Y^-\^<  lOOmV)  at  very  small 
longitudinal  electric  fields  in  the  channel.  The  channel  inversion  charge  density  can 
either  be  measured  using  the  approximation 

QiNV  =  Cox(Vg  -  Vgt)  (3.5) 

or  by  two  ac  capacitance  techniques  to  be  described.  Although  the  value  of  C^,^  in  (3.5) 
can  be  accurately  determined  using  large  area  capacitors,  the  threshold  voltage,  V^j,  is 
rather  ill-defined  in  weak  inversion  near  Yqj,  making  (3.5)  very  inaccurate. 
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A  more  accurate  method  for  extracting  Qjnv  known  as  the  'Split-CV  method, 
was  introduced  by  Koomen  [11].  He  used  it  only  to  study  interface  states  in  the 
inversion  range.  In  this  method  the  gate  capacitance  contributions  of  the 
source/drain/inversion  layer  (Cj+Cj)  are  separated  ('split')  from  those  of  the  bulk.  C^, 
shown  the  in  Figure  3.2  and  they  can  be  easily  measured  using  the  standard  3-terminal 
capacitance  meter.  The  assumptions  are  that  the  small  signal  excess  majority  carrier 
charge  is  supplied  by  the  bulk  current  dig,  whereas  the  contributions  from  interface 
state,  Auj-p,  and  excess  minority  carrier  inversion  charge.  An.  are  supplied  by  dl^  and 
dig.  Thus, 

C,  +  C,  =(l/a))(dlD  +  dlsVdVG  (3.6) 

Cb  =  (l/co)dVdVG  (3.7) 

If  the  interface  charge  density  is  low,  then  C^+Cj  maybe  integrated  over  Vq  to  give 
Qif^  from  weak  to  strong  inversion, 

kc,  +  CJ  dVc  =  qN,Nv  =  QiNV  (3.8) 

This  Koomen  method  was  first  implemented  by  Sodini,  Ekstedt  and  Moll  [52]  to 
measure  the  effective  conductivity  or  linear  mobility,  fj.,i„. 

The  exact  method  for  measuring  the  conductivity  mobility  in  a  surface  inversion 
channel  was  suggested  and  demonstrated  by  Shiue  and  Sah  [53].  They  employed  the 
frequency  dependencies  of  the  ac  conductance  and  capacitance  of  the  channel  treated 
as  a  distributed  one-dimensional  transmission  line.  They  demonstrated  consistency  in 
the  linear  or  effective  conductivity  mobilities  calculated  from  the  experimental  input 
admittance  data  of  the  MOS  channel  transmission  line  in  the  subthreshold  and 
inversion  ranges.  This  exact  method  involved  measuring  the  input  admittance 
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Figure  3.2     Split  -  CV  or  low  frequency  C^^^-Vq  measurement  setup  [33]. 
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(capacitance  and  conductance)  of  the  surface  channel  and  analyzing  the  data  using  a 
distributed-transmission  line  model  known  as  the  circuit  technique  for  semiconductor 
analysis  (CTSA)  [54.  pp.  160-164].  The  CTSA  equivalent  circuit  for  an  MOST,  Fig. 
3.3a  [54,  p.  164].  is  synthesized  using  a  small-signal  expansion  of  the  Shockley 
equations.  The  Cd,_„,,,,i^p  and  C,g.„,.„,^  are  the  drain  and  source  overlap  capacitances. 
Cde  (y)  is  the  distributed  capacitance  of  the  depleted  layer  of  the  surface  space-charge 
layer.  The  charge  storage  capacitance  of  the  electrons  in  the  inversion  channel,  C„(y), 
is  related  to  the  electron  concentration  (number  per  unit  volume).  N(x),  by 

C,-J(q'/kBT)N(x)dx  =  (q^/kBT)NiNv    (F/cm^)  (3.9) 

where 

Njf^^  =  JN(x)dx      (#/cm2)  (3.10) 

Similarly,  the  conductance  element  of  the  channel  is  defined  by 

G^  =  Jq[i (x)N(x)dx  (Siemen  or  mho)  (3.11) 

=  l^Ua'qNlNV 

^i^„  =  I^(x)N(x)dx/NINv  (3.12) 

which  provides  the  rigorous  proof  of  the  definition  given  by  (3.3). 

To  find  C^  and  hence  Nu^y,  the  CTSA  circuit  was  simplified  by  ac  short- 
circuiting  the  gate,  source  and  substrate  (Fig.  3.3b).  The  distributed  capacitance 
element  C(y)  and  the  transmission  line  equations  relating  the  small-signal  input  (v^,  i^) 
to  the  output  (Vj,,  i^)  are, 

C  =  C,(C,,p+C„,)/[C„+(C,,p+CJ].  (3.13) 

Vj  =  VjCOsh(')l.)  -I-  Z^\^svOi]\(C{L),  (3.14) 

id  =  (vyZ„)sinh(YL)  +  i,cosh(il.),  (3.15) 
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Figure  3.3a  and  3.3b     Distributed  Transmission  Line  or  CTSA  model  [32]. 

(a)  Gate  and  Substrate  ac  short  -  circuited  and 

(b)  Gate,  Source  and  Substrate  ac  short  -  circuited 
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Z,  =  [l/jcoCZG„Z]i^  =  (q^/kBTW)N2^^,„,  (3.16) 

Y=  [jcoC/Gji'2  ^  [jco(q/kBT)/^,j''2  =  aw/D,,„)^'2  (3^7) 

where  W  is  the  channel  width,  L  is  the  channel  length,  Z^  is  the  characteristic 
impedance,  y  is  the  propagation  constant,  D,-^=(kgT/q)|i,ij,  is  the  linear  or  effective 
diffusivity,  and  CO  is  the  measurement  frequency.  The  measured  short  circuit 
capacitance,  C^^,,  approaches  a  value  of  CLW/3  when  the  measurement  frequency  is 
much  smaller  than  the  characteristic  frequency  given  by, 

f^  =  (GW/L)/37t(CLW/3)  =  (l/7rL-)(GyC)  -^  Di^^/CttL^).  (3.18) 

where  the  asymptotic  expression  is  for  weak  inversion  when  C^  «  Cj  +  C^^.  If  the 
depletion  and  oxide  capacitances.  Cj  and  C^^,  can  be  obtained  from  the  CV  curve  of  a 
large  area  capacitor,  then,  the  value  of  C„  and  Njj^  can  be  determined. 

It  is  noted  that  by  dc  grounding  the  bulk  of  the  transmission  line  and  ac  short- 
circuiting  the  source  and  drain,  the  input  capacitance  between  the  gate  and  source/drain 
is  identical  to  Cj+C^  in  Koomen's  analysis.  This  provides  a  rigorous  theoretical  proof 
for  Koomen's  empirical  assumptions. 

A  typical  split-CV  or  Cj^^-Vq  measurement,  between  gate  and  the  source/drain 
(base  and  emitter  are  dc  grounded)  of  an  n-channel  BiMOST  device,  is  given  in  Figure 
3.4.  The  integrated  Equation  (3.8)  gives  the  inversion  charge  density  versus  the  gate 
bias  Vq  which  is  also  shown  in  Figure  3.4.  Figures  3.5  and  3.6  illustrate  the 
corresponding  prestress  linear  and  saturation  current-voltage  characteristics  of  the  n- 
channel  BiMOST.  The  drain  current  saturation  condition  is  given  by  Vp  =  V^  >  Vj).sat 
=  V(3  -  Vqj.  By  extracting  the  drain  conductance  from  the  linear  Id-Vq^  for  three 
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Figure  3.4     Initial  prestress  Split-CV  measurement  @10kHz  for  an 
n-channel  BiMOST  with  gate  area  =  1.0E-4cm2  and  an 
oxide  thickness  Xox=150A    Temperature  =  295K. 
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Figure  3.5     Linear  and  logarithmic  plot  of  the  Iq-Vg  and  g  m-^G  characteristics 
for  an  n-channel  BiMOST  with  gate  area  =  1  .OE-4  cm^  and  an 
oxide  thickness  Xqx  =  150A  at  Temperature  =  295K. 
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Sah-Pao  [35]  model  to  include  bulk  charge  effects. 
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values  of  Vps  (90,  100,  1 10m V)  the  prestress  effective  conductivity  or  linear  mobility 
(Fig.  3.7)  can  be  calculated  using  equation  3.2. 

The  second  effective  mobility  investigated  in  this  thesis  was  the  mobility  in  the 
drain-current  saturation  range,  [i^^,.  It  was  extracted  from  the  experimental  I^.^^,  vs 
^G-sat  data  by  least-squares-fit  (LSF)  the  data  to  the  Sah-Pao  MOST  d.c.  equation  with 
bulk  charge  [55].  The  two  LSF  parameters  were  li^^,  and  Wqj  and  the  Sah-Pao 
equations  are: 

iD-sat  =  (li.atCo(W/L){(Vy2)  +  VbVds[  l+(VDs/2(t)p)]^'^ 

-  (4(t)FVB/3)[(l+VDs/2(l)p)2''  -  1]}  (3.19) 

Vds  =  VG_.at-VGT+VB+(Vi/4(t)p)-VB[l+(V^/16(t)^)+(Vc-VGT+VB)/2(^p]i'2  (3.20) 

V3  =  [4qN^e,e„(l)p]"2/C„,  (3.21) 

<1)f  =  kBTlog,(N^/ni)  (3.22) 

where  N^^  is  the  acceptor  concentration  in  the  silicon  and  r^  is  the  intrinsic  carrier 

concentration.  The  temperature  dependence  of  n^  is  given  by  [54,  p. 52] 

n,  =  2.50939xl0'^(m*/m)^'2(j/3Q0)3/2gxp[_£^(2kBT)]  (3.23) 
where  the  energy  gap.  E^,  and  the  effective  mass  ratio,  m*/m,  are  given  by 

Eg  =  1.1700  +  (1.059X10-^)T-  (6.05X10-^)T2     (T<150°K)  (3.24) 

Eg  =  1.1786 -(9.025X10-^)T-(3.05X10-'')T2     (T>150°K)  (3.25) 
m*/m  =  (mj^imp/m^)^'^ 

=  0.81577+3.4353X10-'T[l-(T/437.6)+(T/814.2)2+(T/1356)l  (3.26) 

mf^/mp  =  1 .065/0.647  ^  f(T)  (3.27) 
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Figure  3.7     Drain  (or  channel)  conductance,  gj,  and  conductivity 
effective  electron  or  linear  mobility,  JLiun,  as  a  function  of 
the  d.c.  gate  bias  voltage,  Vq. 
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Deionization  of  the  dopant  impurities  at  low  temperatures  is  also  taken  into  account  by 
solving  the  following  quadratic  for  the  majority  carrier  hole  concentration,  P  of  an  n- 
channel  BiMOST, 

p2  +  PKa  -  KN^  =  0  (3.28) 

Kd  -  (l/gA)Nvexp[-EA/kBT]  (3.29) 

Nv  =  2.50X10''^  {mp/my\T/300y-^  (3.30) 

where  Ny  is  the  effective  density  of  states  in  the  valence  band,  and  mp  is  the  density  of 
states  effective  mass  in  the  valence  band.  The  degeneracy  factor,  g^,  is  assumed  equal 
to  4  and  the  activation  energy  for  boron  is  taken  to  be  45meV.  Similar  expressions  are 
used  for  phosphorus  in  the  p-channel  case  with  E^  replaced  by  Ej)=44meV. 

The  shift  in  the  gate  threshold  voltage,  AVq-j-,  as  extracted  from  the 
aforementioned  two-parameter  fit  to  the  drain-saturation-current  data  will  be  plotted 
unless  otherwise  noted  since  V^-p  can  also  be  extracted  from  the  linear  I^-Vq 
characteristics. 

As  indicated  above,  fi]j„  was  obtained  at  Vp^  — >  0  and  Ey  ^  0  so  there  was  no 
hot  electron  effect  or  mobility  reduction  with  increasing  longitudinal  electric  field. 
However,  [l^^^  was  obtained  at  Vq  =  Vq  =  power  supply  voltage,  hence  hot  electron 
effects  may  be  important  in  very  short  channels.  This  was  not  included  in  the  Pao-Sah 
equation  (3.19),  but  can  be  readily  modified  [10,  pp. 667-670]  for  short  channels. 


CHAPTER  4 
OXIDE  TRAP  CHARGING  AND  DISCHARGING  EXPERIMENTS 


Chapter  3  reviewed  the  fabrication  and  operation  of  the  very  versatile  BiMOST 
test  structure  and  mobility  extraction  methods  for  the  two  effective  mobilities,  ^ijj^  and 
fij^j.  In  this  chapter  the  effects  of  charging  and  discharging  traps  in  the  gate  oxide  layer 
on  the  two  effective  mobilities,  the  drain  saturation  current,  and  the  threshold  voltage 
will  be  demonstrated  by  experiments.  The  injection  mechanism  and  the  nature  of  the 
traps  being  charged  and  discharged  during  the  stress  cycles  will  also  be  discussed.  The 
results  in  this  chapter  are  obtained  on  the  n-channel  BiMOST  and  are  therefore  related 
to  electron  conduction  in  the  surface  channel.  The  results  support  the  theory  of 
increased  Coulombic  scattering  as  the  net  trapped  oxide  charge  density  increases,  and 
decreased  scattering  as  the  net  trapped  oxide  charge  density  decreases.  Furthermore, 
the  changes  in  M-un,  M-saf  lo-saf  and  V^t  ^^^  shown  to  be  consistent  with  the  expected 
sign  of  the  charging  and  discharging  mechanism.  The  linear  proportionality  of  the 
mobility  change  with  positive  and  negative  oxide  charge  density  suggests  a  single 
species  of  oxide  trap  with  three  charge  states  (+,  0,  -)  as  anticipated  from  the  bridging 
oxygen  vacancy  or  E'  center. 

Since  the  I^.^^^,  measurement  is  taken  at  the  condition  of  V^  =  V^  >  V^  -  V^j,  the 
MOST  may  actually  suffer  CHEi  stress  during  the  poststress  characterization.  To  limit 
any  such  measurement  stress,  the  measurement  voltages  were  kept  low  and  only  a 
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small  number  of  triplet  data  points,  with  a  minimum  number  of  samples  for  noise- 
suppressing  averaging,  were  measured  at  each  of  the  predetermined  voltages  shown  in 
Figure  3.6.  To  insure  that  the  device  is  not  inadvertently  stressed  during 
measurements,  the  linear-  and  saturation-IV  and  the  split-CV  were  measured 
repetitively  at  the  intended  measurement  voltages  on  an  unstressed  MOST  to  detect 
any  shift  in  Wqj  during  future  poststress  measurements.  The  results  are  shown  in 
Figure  4. 1  where  the  shift  in  the  gate  voltage  was  measured  at  a  constant  drain  current 
of  l|iA  and  plotted  versus  the  measurement  cycle  to  a  total  of  72  cycles.  In  the  first  22 
cycles  only  the  linear  IV  and  split-CV  were  measured  while  the  remaining  50  cycles 
included  the  saturation  IV  measurement.  Figure  4.1  shows  that  AV^j  is  less  than 
0.2mV  or  the  oxide  charge  changes  by  less  than  3x10^  q/cm^  (from  AQq^  =  -  C^AV^x 
where  C„  =  Eq/^ox  i^  the  oxide  capacitance)  after  72  cycles  at  an  experimental 
resolution  better  than  ±50|LiV  (±10*  q/cm^).  Thus,  this  test  assures  accuracy  of  the 
degradation  data  since  the  test  shows  that  the  MOST  was  not  stressed  significantly 
during  the  measurement  part  of  the  stress-and-measure  (SAM)  experiments  on 
nBiMOSTs  and  nMOSTs  undertaken  in  this  thesis  and  described  in  the  following 
sections  and  the  next  chapter. 

4. 1  Net  Negative  Charge  Build-Up  Due  to  SHEi 

The  first  stress-and-measure  or  SAM  experiment  was  to  stress  the  BiMOST  by 
SHEi  at  a  constant  5MV/cm  gate  oxide  electric  field.  The  gate  current  was  kept 
constant  during  the  stress  at  lG=3nA  and  the  reverse  bias  on  the  collector  base  junction 
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Figure  4. 1     Ths  shift  in  the  gate  voltage,  AV^,  of  an  n-channel  BiMOST  due 
to  repetetive  measurement  of  Id-Vq,  Split  -  CV,  and  saturation. 
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was  kept  constant  at  4V  (Vcb. breakdown  >  15V).  The  shift  of  the  following  MOST 
characteristics  and  material  parameters  due  to  the  SHEi  stress  were  measured:  linear 
and  saturation  Iq-Vg,  transconductance  g^^,  split-CV,  Qjj^,  |i,,„  and  [i^^^,  and  g^.  They 
are  plotted  as  a  function  of  V^  in  Figures  4.2a-d  at  two  stress  fluences,  Nj^j  =  0.  1.3 
and  2.7x1 0'^cm"^.  The  injection  stress  fluence  is  defined  as  the  integrated  electron 
current  flowing  through  the  gate  oxide, 

Nij;j  =  (l/qA)IlGdt      (electron/cm^-sec)  (4.1) 

where  A  is  the  gate  area  and  1^  is  the  d.c.  gate  current.  In  Figure  4.3a  the  effective 
conductivity  or  linear  mobility  is  replotted  as  a  function  of  the  inversion  charge  density 
QiNv  which  is  proportional  to  the  effective  field  in  the  silicon  surface,  as  derived  in 
chapter  2.  thus,  this  figure  illustrates  the  true  mobility  degradation  by  eliminating  the 
effect  of  AVgj.  For  comparison,  the  approximate  linear  mobility  is  also  computed 
from  the  approximated  inversion  charge  density  using  (3.5).  Q^j^y  =  C^^CVq  -  V^j). 
The  approximate  [i^^  versus  the  approximate  Qjj,^  is  plotted  in  Figure  4.3b.  The  value 
of  Yqj  used  here  was  extracted  from  the  linear  Id-Vq  i^^ar  the  peak  g^  which  is  the 
industry  standard  criteria.  Comparing  Figures  4.3a  and  4.3b.  it  is  evident  that  the 
approximate  Qj^v  underestimates  ii^^  due  to  overestimating  Qu^y  from  using  (3.5), 
nevertheless,  the  decreasing  trend  in  the  mobility  with  increasing  stress  fluence.  Njj^tj,  is 

G 

Still  clearly  evident  in  Figure  4.3b.  These  results  are  representative  of  the  X^j^^lVOA 
BiMOST  with  a  gate  area  of  A=100xlOO|-im^.  The  positive  gate  voltage  shift  indicates 
a  buildup  of  net  negative  charge  trapped  in  the  oxide.  The  extracted  threshold  voltage 
shift,  AVqj,  is  shown  in  Figure  4.4  versus  the  electron  injection  fluence  through  the 
gate.  The  repeatability  of  these  results  is  monitored  versus  fluence  rather  than  stress 
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Figure  4.2  a     Shift  in  the  linear  IV  of  an  n-channel  BiMOST  due  to 

SHEi  stress  at  a  constant  Eqx  =  4MV/cm.    The  curves  for 
three  different  fluences  are  plotted. 
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Figure  4.2  b     Shift  in  the  saturation  IV  of  an  n-channel  BiMOST  due  to 
SHEi  stress  at  a  constant  Eqx  =  4MV/cm.    The  curves  for 
three  different  tluences  are  plotted. 
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Figure  4.2  c     Shift  in  the  split-CV  and  Qjnv  of  an  n-channel  BiMOST 
due  to  SHEi  stress  at  a  constant  Eqx  =  4MV/cm.    The 
curves  for  three  different  fluences  are  plotted. 
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Figure  4.2  d     Shift  in  the  effective  mobility  and  gj  of  an  n-channel  BiMOST 
due  to  SHEi  stress  at  a  constant  Eqx  =  5MV/cm.    The  curves  for 
three  different  fluences  are  plotted. 
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Figure  4.3a      Shift  in  the  effective  mobility  of  an  n-channel  BiMOST  due 
to  SHEi  stress  at  a  constant  Eqx  =  4MV/cm.    The  curves  for 
three  different  fluences  are  plotted. 
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Figure  4.3  b  Shift  in  the  effective  mobility  of  an  n-channel  BiMOST  due 
to  SHEi  stress  at  a  constant  Eqx  =  4MV/cm.  The  curves  for 
three  different  fluences  are  plotted  versus  an  approximated  Qinv- 
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Figure  4.4     Threshold  voltage  shift  of  n-channel  BiMOST  during  SHEi  at 
T=295°K,  Eox=  4MV/cm,  Vcb  =  4V  and  at  a  constant  lo^SnA. 
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time  since  the  gate  current  sometimes  varies  from  device  to  device.  In  Figure  4.5a  and 
4.5b  the  degradation  of  the  saturation  current.  AIp.,^,.  is  plotted  versus  the  injection 
fluence,  Nj-^j.  In  Figure  4.5a  the  value  of  Ij)_^^i  is  extracted  at  a  constant  value  of  Wq  = 
^G  >  ^D-sat  -  ^G  ~  ^GT  whereas  in  Figure  4.5b  I^.^^,  was  extracted  at  V^  -  V^x:  which 
includes  the  variation  of  V^-j.  with  NI^;J.  Thus.  Figure  4.5a  illustrates  the  contributions 
to  AIjj.j^j  from  both  the  shift  in  the  threshold  voltage,  AV^-p,  and  the  degradation  or 
change  of  the  saturation  mobility  Afi^^,  whereas  Figure  4.5b  AId,,^,  contains  only  the 
contributions  of  Afi,^,.  The  percent  change  in  the  effective  saturation  mobility,  [i,^,.  and 
the  effective  conductivity  or  linear  mobility,  [ly^^,  versus  the  injection  tluence  is  plotted 
in  Figure  4.6.  The  value  of  \i^j^  was  monitored  for  a  variety  of  constant  Qjvjy  values  to 
maintain  a  constant  surface  field  in  the  silicon.  As  clearly  evident  in  Figures  4.3a  and 
4.3b,  the  mobility  curves  converge  at  large  Q^^y  and  hence  the  shift  in  fi,,„  is  less  as 
shown  in  Figure  4.6.  In  Figure  4.6.  the  prestress  values  of  liji^  range  from  557cmVV- 
sec  near  the  peak  of  the  curve  at  QINv=1.34xI0'^C/cm^=8.4xl0"q/cm^  to  532cm W- 
sec  atQjf^=2.51xlO"^C/cm^=1.57xlO'^q/cm^,  while  the  prestress  value  of  the  effective 
saturation  mobility  is  518cmVV-sec.  Figure  4.6  indicates  that  the  saturation  mobility 
is  more  sensitive  to  the  increasing  charge  in  the  oxide.  To  summarize  the  results, 
Figure  4.7  plots  the  extracted  AVqj,  A^ij-^,  Aji^g,  versus  the  total  AI^.^^,  (I^.s^j  extracted 
Vp  =  Vg  =  IV  =  constant). 

In  order  to  quantify  the  effects  of  interface  traps  on  the  mobility,  in  addition  to 
the  just-described  effects  from  the  oxide  charge  buildup  during  the  stress,  two  methods 
were  employed  to  measure  the  increase  in  the  interface  trap's  quantum  density-of-state, 
Dj-p.  The  first  is  the  commonly  used  method  of  monitoring  the  change  of  the 
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Figure  4.5  a     Shift  in  the  saturation  current  of  an  n-channel  BiMOST  due 
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Figure  4.5  b     Shift  in  the  saturation  current  of  an  n-channel  BiMOST  due 

to  SHEi  stress  at  a  constant  Eox=  4MV/cm.    The  value  of  Ip.sat 
was  extracted  for  a  varying    value  of  ^D-sai  =  ^G-sat  plus    the 
extracted  value  of  Yqj. 
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Figure  4.6       Percentage  shift  in  the  conductivity  effective  mobility,  [i]^,,, 
for  varying  values  of  Q^f^y,  and  in  the  saturation  effective 
mobility,  |isa,  plotted  versus  fluence  for  an  n-channel  BiMOST 
during  SHEi  stress  at  constant  Eox=4MV/cm  and  Vcb  =  4V. 
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Figure  4.7      The  relative  change  in  threshold  voltage,  Vq-j-,  the  conductivity 

effective  mobility,  |i]i„  and  the  saturation  effective  mobility,  |J.sat, 
versus  the  corresponding  degradation  in  lD_sat  for  an  n-channel. 
BiMOST  stressed  at  Eqx  =  4MV/cm. 
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subthreshold  voltage  swing,  AS,  in  which  the  inverse  slope  in  the  prestress 
subthreshold  portion  of  the  Id-V^s  data  is  measured  and  compared  to  that  measured 
after  each  stress.  Although  the  ideal  theoretical  value  at  room  temperature  is 
AVc=2.303kgT/q=60mV  per  decade  of  I^,  a  measured  value  of  80mV  per  decade  is 
more  common  due  to  thick  oxide  and  high  substrate  doping  variation  [10,  Fig. 682.1  on 
p. 655].  The  relation  of  the  density  of  the  generated  interface  trap  to  the  subthreshold 
swing,  AS,  is  given  by  {See  reference  [10],  Eqns.(683.10)  and  (683.11)  on  p.666.} 

ADiT  =  D,T(t,,,ess)  -  Dit(O)  =  (Cox/q)(q/2.303kBT)AS  (tt/eV-cm^).  (4.2) 

In  our  samples,  the  prestress  density  of  interface  traps  is  negligible  and  may  be 
assumed  zero, 

AS  =  S(Dit)  -  S(Dix=0)  =  S(Dit).  (4.3) 

The  second  approach  is  a  new  method  for  measuring  Djj  buildup  recently 
proposed  by  Neugroschel  and  Sah,  and  demonstrated  by  them  and  their  graduate 
students  [12].  This  method  takes  advantage  of  the  unique  feature  of  the  vertical  bipolar 
junction  transistor  (BJT)  structure  in  the  BiMOST  and  the  fact  that  the  buildup  of  the 
BIT  base  recombination  current,  Alg,  is  directly  and  solely  proportional  to  the  buildup 
of  the  interface  trap  density  because  the  additional  recombination  of  the  minority 
carriers  (electrons  in  the  p-base  layer  of  n/p/n  BJT  of  the  nBiMOST)  occurs  at 
generated  interface  traps.  They  demonstrated  the  extraction  of  the  interface  trap 
density  from  the  change  of  the  stress-induced  increase  in  the  base  current  after  each 
stress  cycle.  This  new  method  was  coined  and  acronymed  by  Sah  as  the  so-called 
direct-current  current-voltage  measurement  or  DCIV. 
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The  relationship  between  the  change  in  the  base  current,  Alg,  the  surface 
recombination  velocity.  S„,  and  the  emitter  base  bias,  Vgg,  is  given  by  [12] 

AIb  =  (qAn,ASo/2)exp(qVBE/kBT)  (4.4) 

ASo  =  (7t/2)ao0,hANn-  (4.5) 

where  Gq  is  the  cross-section  of  carrier  capture  at  the  interface  traps,  0;^  is  the  carrier 
thermal  velocity,  and  ANj-j-  is  the  stress-generated  interface  trap  density  (#/cm^).  The 
base  current,  Ig,  versus  the  gate  to  base  d.c.  bias,  Vcg.  measured  during  the  SHEi  stress 
at  Eox=5MV/cm  is  shown  in  Figure  4.8  for  three  stress  fluences,  Nj^j  =  0,  1.3  and 
2.7xl0*^cm"^.  A  comparison  of  the  two  measurements,  AS  and  %AIb,  is  shown  in 
Figure  4.9  versus  the  injection  fluence  where  %AIB/100  is  calculated  from, 

%Ay  100  =  [lB(t,ttess)  -  Ib(0)]/Ib(0)  (4.6) 

where  Ib(0)  is  the  prestress  value  and  %^^^^  is  the  stress  time  or  the  stress  fluence. 

The  results  of  the  two  interface  trap  monitoring  methods  indeed  track  each  other 
quite  well.  The  value  for  Dj-j-  plotted  was  calculated  from  AS  assuming  there  are  no 
interface  traps  in  the  unstressed  sample.  A  calculation  of  Dj^  from  Alg  would  require 
an  assumption  for  the  cross-section,  Gq.  In  the  following  chapter  a  more  thorough 
comparison  of  AS  with  AIb  will  be  given  to  separate  the  effects  of  interface  traps  and 
areal  (y-direction)  nonuniformity  of  trapped  oxide  charge. 

In  order  to  investigate  the  effect  of  different  charge  distributions  within  the  oxide 
(x-direction)  another  device  from  the  same  170 A  group  was  stressed  at  a  lower  gate 
oxide  electric  field  of  4MV/cm  with  all  other  stress  and  characterization  parameters 
remaining  the  same.  Minority  carriers  traveling  through  the  gate  oxide  during  the 
SHEi  will  now  have  a  lower  kinetic  energies  and  will  give  a  different  oxide  charge 
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Figure  4.8        Increase  of  the  base  current    of  an  n-channel  BiMOST  due 
to  the  build-up  of  interface  states  during  SHEi  stress  at  a 
constant  Eox=4MV/cm.  Ib  is  plotted  for  three  different  fluences. 
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Eox=4MV/cm. 
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distribution  in  the  thickness  direction  or  x-direction  of  the  thin  gate  oxide  film,  Qqj  = 
Qot(x).  Furthermore,  the  lower  electric  oxide  field  of  4MV/cm  will  give  a  lower  field- 
aided  and  tunneling  emission  rates  of  electrons  trapped  in  shallow  energy  levels  in  the 
oxide.  Figures  4.10-4.14  give  an  illustration  of  the  oxide  electric  dependence  of  the 
stress  induced  AVqj,  AIq.,^,,  Ajiji^,,  Ap,,^,,  AS  and  Alg  versus  fluence  for  Eox=3MV/cm 
and  Eox=4MV/cm.  The  degradation  of  all  these  parameters  indeed  increases  with 
increasing  gate  oxide  electric  field.  Furthermore  the  lower  field  stress  appears  to  cause 
a  saturation  of  these  parameters  at  a  lower  fluence,  indicating  that  the  oxide  trap 
charging  or  discharging  rate  may  be  lower  and  that  the  traps  being  filled  by  or  emptied 
of  electrons  have  a  larger  capture  or  emission  cross-section.  In  order  to  make  a  fair 
comparison  of  the  degradation  in  the  extracted  AHy^  and  Afi^^^  at  the  two  gate  oxide 
electric  fields.  Figure  4.12  should  be  replotted  versus  the  trapped  charge  density,  Qqt- 
Hence  if  the  spatial  distribution  (in  x-direction  or  the  oxide-thickness  direction)  of  the 
stress  built-up  oxide  traps  at  different  Eqx  and  stress  conditions  were  the  same,  the 
degradation  in  the  two  effective  mobilities  would  be  the  same.  However,  the  two 
mobilities  and  drain  current  are  also  degraded  and  the  gate  threshold  voltage  is  also 
changed  by  the  stress-generated  interface  traps  which  carmot  be  reversed  during  the 
electronic  neutralization  of  the  oxide  traps  by  electron  thermal  capture  and  tunnel 
emission,  because  the  generated  interface  traps  cannot  be  removed  or  annealed 
electronically.  The  two  experiments  described  in  the  next  two  subsections  will 
demonstrate  the  reversibility  of  !!„„,  [i,^,,  lo-saf  and  Y^t  during  the  positive  and  negative 
charging  and  neutralization  cycles  of  the  oxide  traps. 
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Figure  4.10     Threshold  voltage  shift  of  n-channel  BiMOST  during  SHEi  at 
Eox=3MV/cm  and  Eox=4MV/cm;  Vcb=4V;  T=295°K. 
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Figure  4.11a    Shift  in  the  saturation  current  of  an  n-channel  BiMOST  due 

to  SHEi  stress  at  Eox=3MV/cm  and  Eox=4MV/cm.    The  value 
of  lo.sat  was  extracted  at  a  constant  value  of  Vi5-sat=  ^G-sm- 
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Figure  4.11  b    Shift  in  the  saturation  current  of  an  n-channel  BiMOST  due 

to  SHEi  stress  at  Eox=3MV/cm  and  Eox=4MV/cm.  The  value 
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Figure  4.12     Percentage  shift  in  the  conductivity  effective  mobiUty,  [i^^,  and 

saturation  effective  mobility,  [i^^i  of  electrons  plotted  versus  fluence 
for  an  n-channel  BiMOST  during  SHEi  stress  at  Eox=3MV/cm  and 
Eox=4MV/cm  and  Vcb=4V. 
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Figure  4.13      Shift  in  the  subthreshold  slope,  S,  and  the  corresponding  increase 
in  the  base  cun-ent,  AIb,  of  an  n-shannel  BiMOST    due  to  the 
the  build-up  of  interface  states  during  SHEi  stress  at  a  constant 
Eox=3MV/cm  and  Eox=4MV/cm  and  Vcb=4V. 
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Figure  4.14     Percentage  shift  in  the  conductivity  effective  mobility,  ii^eff  and 
saturation  effective  mobility,  \i^^i  of  electrons  plotted  versus  the 
net  trapped  charge  density  for  an  n-channel  BiMOST  during  SHEi 
stress  at  Eox=3MV/cm  and  Eox=4MV/cm;  Vcb=4V. 
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4.2  Charging  and  discharging  the  deep  oxygen  vacancy  center  at  Er  -  7eV 

One  of  the  recently  identified  and  studied  defects  in  the  SiOj  quartz,  optical 
fibers,  and  MOST's  gate  oxide  is  the  E'  center,  commonly  known  as  the  bridging 
oxygen  vacancy  center,  pictured  in  Figure  4.15.  The  chemical  symbol  to  be  used  in 
this  thesis  will  be  Vq  for  the  bridging  oxygen  vacancy  center,  and  Vq,  Vq,  and  V3  for 
the  center  in  the  three  charged  states  as  proposed  by  Sah  to  account  for  the  electronic 
transitions  between  the  multi-charge-states  of  the  oxygen  vacancy  center  [54.  Fig.B3.4 
on  page  422.].  Thompson  [56,  57]  monitored  the  buildup  kinetics  of  positive  charge  in 
the  oxide  during  FNTE  or  SHE  injection  stress  at  high  oxide  electric  fields  (Fig. 
4.15c).  His  findings  were  attributed  by  Sah  to  the  bridging  oxygen  vacancy  [54, 
Fig.B3.4  on  page  422.].  The  E'  center  has  been  well  known  in  fiber  and  optical  glass 
investigations  and  its  microscopic  configuration  indicated  in  Figure  4.15d  was 
delineated  by  electron  paramagnetic  resonance  spectroscopy  [58].  It  is  but  only  one  of 
the  charge  state  configurations  of  the  oxygen  vacancy  center,  that  is  detectable  by  spin 
resonance  of  an  unpaired  electron.  The  multi-charge-state  and  transition  energy  band 
model  proposed  by  Sah,  as  depicted  in  Figure  4.15c,  suggests  that  positive  oxide 
charge  is  created  when  one  of  the  bound  electrons  at  an  initially  neutral  oxygen 
vacancy  site  [(Si-0)3=Si-Si=(0-Si)3],  is  impact  emitted  or  released  by  an  energetic 
electron  (kinetic  energy  >  7eV).  In  Thompson's  experiment  suggested  by  Sah, 
electrons  were  injected  into  the  gate  oxide  from  the  gate  electrode  during  FNTEi  or 
from  the  substrate  during  SHEi.  They  were  then  accelerated  to  7eV  by  the  oxide 
electric  field  to  gain  sufficient  kinetic  energy  to  impact  release  the  bond  electron  at  the 
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Figure  4.15     Energy  band  diagram  illustrating  the  various  states  of  the  E'  center 
otherwise  known  as  the  oxygen  vacancy;  (a)  thermal  or  tunneling 
emission  from  a  shallow  state,  (b)  thermal  capture  at  a  shallow  state, 
(c)  generation  of  positive  charge  by  impact  emission  of  a  bound 
electron  by  an  energetic  electron  injected  from  the  substrate  or  poly- 
gate,  and  (d)  theiTnal  capture  at  the  deeper  level  (Ec-Et=7eV) 
[51,Fig.  B3.4] 
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neutral  oxygen  vacancy.  These  positively  charged  oxygen  vacancy  centers  can  then  be 
neutralized  by  capturing  electrons  (Fig.  4. 15d)  from  low-field  SHEi.  at  sufficiently  low 
oxide  electric  field  such  that  the  kinetic  energy  gained  by  injected  electrons  in  the 
oxide  is  less  than  about  7eV  and  impact  emission  of  the  newly  captured  electrons 
would  be  avoided. 

The  results  of  charging  and  discharging  the  7eV  oxygen  vacancy  center  under 
high  and  low  oxide  electric  fields  and  the  effect  of  charge  on  the  electron  mobility  and 
^D-sat  '^ill  °ow  ^^  discussed.  An  n-channel  BiMOST  (Xqx=170A)  was  first  subjected 
to  the  low  field  EQx=1.5MV/cm  SHEi  stress  to  fill  any  already  positively  charged 
states  in  the  oxide.  The  positive  threshold  voltage  shift.  AVq^,  is  shown  in  Figure 
4.16a  and  enlarged  in  Figure  4.16b  at  low  fluence  (Njjjj<4X10''*/cm^).  Next  electrons 
were  Fowler-Nordheim  tunnel  injected  from  the  poly-gate  into  the  oxide  to  impact  emit 
the  trapped  electrons  from  the  E^.  -  7eV  level  of  the  Vq  center  in  the  oxide  giving  Vq  -i- 
e"*  ^  Vq  -I-  2e".  The  resulting  negative  threshold  voltage  shift  is  illustrated  in  Figure 
4.16a.  Finally  the  oxide  was  again  SHE  injected  at  EQx=l-5MV/cm  to  refill  the 
positively  charged  centers,  Vq  -t-  e"  ^  Vq,  as  indicated  in  the  enlargement  of  Figure 
4.16c  for  (9x10^^  <  Njj^j  <  9.4x1 0'^^cm"^).  Clearly  V^  recovers  almost  completely,  to 
within  lOmV  of  the  initial  value  after  low-field  SHEi. 

In  Figures  4.17  and  4.18  the  corresponding  Al^.s^t  are  shown.  The  value  of  I^.^g, 
was  extracted  at  a  constant  Vd=Vq  in  Figure  4.17  and  at  V^^Vq  adjusted  by  AVqj  in 
Figure  4.18  to  separate  the  mobility  and  threshold  voltage  contributions  to  AI^.j^j. 
These  two  figures  differ  drastically  due  to  the  fact  that  the  built-up  oxide  charge  does 
not  seem  to  affect  the  value  of  \l^^^  as  much  as  V^j.    Therefore  the  degradation 
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Figure  4.16  (a)  Threshold  voltage  shift  of  n-channel  BiMOST  during  SHEi  at 
Eox=  1.5MV/cm  enlarged  in  (b)  for  Ninj<4  x10^\  followed  by  FNTEi 
from  the  gate  at  Eqx  =  8MV/cm,  for  4  xlO^''<  Ninj<  9  xlO'^  ,and  again 
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Figure  4.17     Shift  in  the  saturation  current  of  an  n-channel  BiMOST  during 

SHEi  at  Eox=l-5MV/cm  enlarged  in  (b)  for  Ninj<4  xlO'"*,  followed 
by  FNTEi  from  the  poly-gate  at  Eox=8MV/cm.  for  4  xlO<Njv,j<9  xlO^^ 
and  again  by  SHEi  at  Eox=l-5MV/cm  for  9  x10^^<Ninj<9  xlO^'  in  (c). 
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Figure  4.18     Shift  in  the  saturation  current  of  an  n-channel  BiMOST  during 

SHEi  at  Eox=l-5MV/cm  enlarged  in  (b)  for  Ninj<4x10^^  followed 
by  FNTEi  from  the  poly-gate  at  Eox=8MV/cm,  for  4  x10<Ninj<9  xlO^^ 
and  again  by  SHEi  at  Eox=l-5MV/cm  for  9  x10^^<Ninj<9  xlO^^  in  (c). 
The  value  of  Ii^.^^t  was  extracted  at  a  constant  VD.sat=VG-sat  plus  AVgt. 
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in  I0.JJ,  appears  negligible  if  the  value  of  Vo=Ve+AVGj  is  used  to  extract  the  current  as 
in  Figure  4.18.  This  figure  essentially  illustrates  Al^.^a,  due  to  the  degradation  in 
mobility.  This  is  also  evident  in  the  plot  of  the  mobility  degradation  in  Figure  4.19. 
The  three  hollow  marker  curves  in  this  plot  give  |i,i„  extracted  at  at  three  constant 
values  of  Qixv-  These  curves  support  my  assumption  that  the  positively  charged 
oxygen  vacancy,  Vq,  are  filled  with  electrons  during  the  low  field  SHEi  to  neutrahze 
the  centers,  decrease  the  Coulombic  scattering,  and  increase  the  mobility,  as  I 
anticipated  in  designing  this  set  of  experiments.  On  the  other  hand  during  the  high 
field  FNTEi  the  electrons  are  impact  emitted  from  the  7eV  level  leaving  behind  a 
positively  charged  center.  This  in  turn  increases  the  Coulombic  scattering  and 
decreases  |a,,^.  The  solid  triangle  curve  in  Figure  4.19  represents  the  saturation 
mobility,  \i^^^,  which  appears  to  be  less  sensitive  to  charging  and  discharging  of  this 
oxide  trap  (oxygen  vacancy)  as  suggested  by  the  difference  between  Figures  4.17  and 
4.18.  Thompson  showed  [56]  that  the  centroid  for  this  positive  charge,  Vq,  is 
approximately  80-90A  from  the  injecting  interface  or  at  about  central  plane  of  the 
170A  oxide  layer  investigated  here.  Nevertheless  ]l^^^  does  degrade  slightly  during  the 
long  FNTEi  impact  emission  cycle  and  does  recover  immediately  following  the  start  of 
the  second  low  field  SHEi  cycle.  By  isolating  the  initial  portion  of  the  high-field 
FNTEi  and  plotting  the  extracted  positive  oxide  charge  density.  +Qot  -  ~(Cox/Q)^^gt' 
against  the  degradation  of  [i^^^  (at  a  constant  Qij^v=l-5xlO~'C/cm^=9.4xlO''q/cm^)  an 
essentially  linear  dependence  on  the  oxide  charge,  Qox^  is  revealed  (Fig.  4.19d). 

Finally  it  is  noted  that  the  generation  of  interface  states  as  measured  by  AS  and 
Alg  was  neghgible  during  the  SHEi  portions  of  the  stress  (Figs.  4.20b  and  4.20c). 
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Figure  4.19     Percent  degradation  of  the  mobility  of  an  n-channel  BiMOST  during 
SHEi  at  Eox=l-5MV/cm  enlarged  in  (b)  for  Ninj<4x10'\  followed 
by  FNTEi  from  the  poly-gate  at  Eox=8MV/cm,  for  4  xlO  <  Njnj  <  9  xlO 
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Figure  4.19  d    Effective  electron  mobility  shift  versus  the  impact-emission  generated 
+QoT'  trapped  in  the  gate  oxide  during  the  high-field  FNTEi. 
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Figure  4.20  (a)  Subtheshold  slope  and  base  current  degradation  of  n-channel  BiMOST 
during  SHEi  at  Eox=  1.5MV/cm  enlarged  in  (b)  for  Nj^j  <  4  xlO'"*,  followed 
by  FNTEi  from  the  gate  at  Eqx  =  8MV/cm,  for  4  xlO'^<  N,nj  <  9  xlO^^  ,and 
again  by  SHEi  at  Eox=1.5MV/cm  for  9  xlO'^<  Njnj  <  9.4  xlO'-^  enlarged  in  (c). 
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However,  some  Dj-j-  are  generated  during  the  FNTEi  portion  of  the  stress. 
Interestingly  though  the  subthreshold  slope  recovers  to  within  10%  of  its  original  value 
during  the  second  SHEi  stage  whereas  the  shift  in  Ig  remains  constant.  This  is  most 
likely  due  to  the  fact  that  the  subthreshold  swing  is  sensitive  to  the  areal  nonuniformity 
of  the  oxide  traps  from  charging  and  discharging  in  addition  to  the  stress-generated 
interface  traps,  whereas  the  peak  value  of  Ig  depends  only  on  the  change  or  increase  of 
the  recombination  rate  from  the  increasing  density  of  the  interface  traps.  Hence  the 
fact  that  AS  recovers  and  Alg  does  not  indicates  that  trap  charging  and  discharging  are 
nearly  areally  uniform  in  these  experiments.  And  the  experiments  demonstrate  again 
that  AS  is  not  as  sensitive  and  unambiguous  as  Alg  in  monitoring  the  interface  traps. 

4.3  Charging  and  discharging  the  shallow  oxygen  vacancy  center  at  Er;  -  leV 

Thompson  [57]  in  his  doctoral  thesis  also  investigated  the  shallow  level 
associated  with  the  oxygen  vacancy  located  at  Eq-  leV  below  the  oxide  conduction 
band.  Since  electrons  trapped  at  this  level  are  easily  thermally  emitted  or  detrapped  at 
room  temperatures  the  leV  level  was  filled  at  77°K  using  low  field  SHEi  stress.  Then 
the  electrons  were  tunneled  out  of  the  trap  to  the  oxide  conduction  band  using  the 
isochronal  electric  field-stimulated  emission  (EFSE)  technique.  In  this  measurement 
the  oxide  electric  field  is  incrementally  increased  by  the  applied  gate  bias  with  the 
source,  drain,  well,  and  substrate  grounded  for  a  predetermined  amount  of  time.  As 
Eqx  increases,  the  trap  to  band  tunneling  probability  increases  until  the  trapped 
electrons  are  emitted  from  the  shallow  traps.  In  this  section  the  effects  of  charging  and 
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discharging  this  shallow  leV  center  on  the  mobility  and  saturation  current  will  be 
illustrated. 

In  this  SAM  experiment,  as  in  the  previous  ones,  it  is  essential  that  the  IV  and 
CV  characterization  measurements  are  limited  to  a  range  of  gate  bias  voltages  lower 
than  that  applied  during  the  stress  in  order  that  the  trapped  electrons  are  not  field 
emitted  from  the  shallow  leV  level  during  the  measurement.  This  cannot  be  attained 
at  room  temperature  because  of  rapid  thermal  emission  of  trapped  electrons  to  the 
oxide  conduction  band  from  the  shallow  level.  Therefore,  in  order  to  obtain  an 
adequately  long  Vq  sweep  for  extracting  |i,j„  while  keeping  the  oxide  electric  field 
sufficiently  low  (<lMV/cm),  an  n-channel  BiMOST  with  a  thick  gate  oxide  was  used 
with  Xox=290A  and  a  circular  geometry  of  W/L=16.8.  After  immersing  the  packaged 
BiMOST  into  liquid  nitrogen  in  the  computer-controlled  dewar,  a  5MV/cm  field  is 
applied  across  the  oxide  for  5  minutes  to  ensure  that  all  the  shallow  leV  levels  are 
empty.  Next  the  BiMOST  is  biased  for  SHEi  at  a  low  gate  oxide  electric  field  of 
IMV/cm.  Since  the  injection  efficiency  is  rather  poor  at  such  a  low  field  the 
collector/base  junction  reverse  bias  was  set  to  7V  and  the  emitter/base  junction  was 
biased  at  3V  initially. 

Figure  4.21  plots  the  shift  in  the  threshold  voltage  versus  the  injection/emission 
time,  since  the  gate  current  is  negligible  during  the  emission  portion  of  the  stress  and 
does  not  contribute  significantly  to  Nj^j.  The  emission  time  for  each  Vq  increment  of 
0.3  V  was  approximately  10  seconds.  In  Figure  4.22  the  emission  portion  of  the  stress 
is  plotted  versus  the  increasing  gate  oxide  field  illustrating  the  Yqj  recovery  as  the 
trapped  electrons  are  emitted.  In  Figures  4.23  and  4.24  the  shift  in  I^^^t  is  plotted  for 
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Figure  4.21  Threshold  voltage  shift  of  an  n-channel  BiMOST  during  SHEi 
stress  at  a  constant  Eox=lMV/cni  and  Vcb=7V,  followed  by 
isochronal  EFSE  for  t>400sec. 
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Figure  4.22     Threshold  voltage  shift  of  an  n-channel  BiMOST  during  isochronal 
EFSE  versus  the  gate  oxide  electric  field  Eox- 
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Figure  4.23      Shift  in  the  saturation  current  of  an  n-channel  BiMOST  due 
to  SHEi  stress  at  a  constant  Eox=  IMV/cm  and  Vcb=7V, 
followed  by  the  isochronal  EFSE  for  t>400sec.    The  value  of 
iD-sat  was  extracted  at  a  constant  value  of  Vo-sat  =  ^o-sat- 
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Figure  4.24      Shift  in  the  saturation  current  of  an  n-channel  BiMOST  due 
to  SHEi  stress  at  a  constant  Eox=  IMV/cm  and  Vcb=7V, 
followed  by  the  isochronal  EFSE  for  t>400sec.    The  value  of 
In-sat  was  extracted  at  a  value  of  Vd.s^^  =  V^.j^,  adjusted  by  AVg^. 
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the  constant  Vd=V(3  and  the  Vd=Vc+AVct  cases  as  before.  Though  the  value  of  Ij),^^^ 
decreases  largely  due  to  the  shift  in  Vgj  there  appears  to  be  a  significantly  larger 
decrease  in  ^1,^^  seen  in  the  Vcx-adjusted  V^ot  (Fig.  4.24).  The  degradation  and  recovery  in 
|i]i„  and  [ij^t  are  shown  in  Figure  4.25  where  |i]i„  is  plotted  for  three  different  values  of 
constant  Qi^v  The  magnitude  of  A|i,^^  is  significantly  smaller  than  that  of  A^i^^,.  This 
may  be  due  to  the  fact  that  the  'Vu-sM=^G-s!it  values  for  extracting  lD..,a,  remain  constant 
during  the  stress  even  though  there  is  a  large  shift  in  V^t  as  shown  in  Figure  4.21. 
This  would  cause  a  significant  change  of  the  voltage  drop  and  longitudinal  electric 
field  in  the  drain  junction  space-charge  region  from  the  carrier  depletion  point  in  the 
channel  to  the  drain  junction,  which  would  alter  the  buildup  rate  of  the  oxide  charge 
and  interface  traps.  Therefore,  the  range  of  current  used  in  the  two-parameter  (Vqj, 
jjj^t)  least-squares-fit  to  the  Sah-Pao  model  is  changed  after  each  stress.  For  example,  a 
large  positive  V^-p  shift  would  have  shifted  the  data  closer  to  the  subthreshold  range 
where  the  Sah-Pao  bulk-charge  model  does  not  apply.  Nevertheless  both  mobilities 
decrease  as  expected  with  the  increasing  trapped  electron  charge  in  the  oxide  due  to  the 
corresponding  increase  in  Coulombic  scattering.  The  mobilities  also  recover  as 
expected  when  the  electrons  are  tunnel  emitted  out  of  the  shallow  traps  leaving  them  in 
the  neutral  charge  state  and  decreasing  the  Coulombic  scattering.  Plotting  the 
degradation  of  the  effective  mobility  during  the  initial  low-field  -Qqt  charging  stage 
versus  the  negative  oxide  charge  density  extracted  from  the  AVq-j.,  a  linear  dependence 
on  -QoT  is  observed  as  shown  in  Fig.  4.26.  This  is  similar  to  that  shown  in  Fig.4.19d 
for  the  +QoT  ^ind  provides  a  second  experimental  proof  of  electron  mobility  reduction 
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Figure  4.25     Percentage  shift  in  the  conductivity  effective  mobility,  fiii^, 
for  varying  values  of  Qintv,  and  in  the  saturation  effective 
mobility,  [ism  plotted  versus  time  for  an  n-channel  BiMOST 
during  SHEi  stress  at  constant  Eox=lMV/cm  and  Vcb  =  7V, 
followed  by  the  isochronal  EFSE  for  t  >  400sec. 


86 


400  I   I    I    I    I    I    I   I    I    I   I    I    I    I    I   I   I    I    I    I    I    I    I    I    I 


> 


> 

< 


300 


200   - 


100   - 


0 


I    I    I    I    I    I    I    I    I    1    I    I    I    I    I    I    I    I    I    I    I    I    I    I    I 


0  -2  -4  -6  -8 

%A|Liii„/(cm'.V^sec-^) 


-10 


Figure  4.26  Shift  of  the  effective  conductivity  electron  mobility  during  the 
low-field  SHEi  filling  of  the  shallow  electron  traps  in  the  gate 
oxide  versus  the  negative  oxide  charge  density. 
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Figure  4.27      Shift  in  the  subthreshold  slope,  S,  of  an  n-channel  BiMOST  due  to 
the  build-up  of  interface  states  and/or  nonuniform  charging  during 
SHEi  stress  at  a  constant  Eox=lMV/cm.and  Vcb=7V,  followed  by 
isochronal  EFSE  for  t  >  400sec. 
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from  Coulomb  scattering  of  the  channel  electrons  by  positive  and  negative  charge 
states  of  oxide  traps  at  the  same  location  in  the  oxide. 

Finally  the  subthreshold  swing  and  the  corresponding  calculated  Djj  are  shown 
in  Figure  4.27.  The  recovery  shows  that  AS  is  most  likely  due  to  areal  nonuniformity 
in  the  oxide  charge  as  discussed  previously  rather  than  true  interface  state  generation 
and  annealing  since  such  a  large  number  (lO'Vcm^)  of  interface  states  would  not  be 
expected  during  low,  IMV/cm.  SHEi  stress,  and  in  addition,  there  is  not  a  viable 
physical  mechanism  to  anneal  out  the  interface  states.  This  is  collaborated  with  the  Alg 
change  during  the  low  field  SHEi  shown  in  Fig.  4.20  which  shows  no  changes  during 
the  "annealing"  cycle  when  AS  reduces  to  zero  while  Alg  remains  high. 

This  data  will  be  expanded  upon  in  the  following  chapter  to  illustrate  the 
separation  of  interface  states  and  areal  charge  nonuniformity. 


« 


CHAPTER  5 

SEPARATION  OF  INTERFACE  TRAPS  AND  AREAL 

OXIDE  CHARGE  NONUNIFORMITY 


The  interface-trap  caused  distortion  of  the  gate  capacitance-voltage  (C^-Vq  or 
CV)  characteristics  of  the  MOSC  [10,  Fig.403.1  on  p.324]  and  the  d.c.  VVq  (IV) 
characteristics  of  the  MOST  in  both  the  subthreshold  and  the  inversion  ranges  [10, 
section  628  on  pp. 653-662]  are  well  documented  in  the  literature.  However  it  is  also 
well-known  that  an  areally  nonuniform  distribution  of  charged  oxide  traps  (Qqt)  will 
distort  the  CV  and  IV  characteristics.  In  addition,  the  FV  characteristic  is  also  distorted 
by  the  channel  mobility  variation  with  Vq  as  discussed  in  chapter  2  and  illustrated  in 
chapter  4.  In  the  last  chapter  it  was  shown  that  the  MOST  subthreshold  swing,  AS,  and 
base  current  change,  Alg,  could  be  used  to  monitor  the  increase  in  interface  traps  but 
that  they  could  also  be  affected  by  areally  (y-direction)  nonuniform  trapped  oxide 
charge.  In  this  chapter  a  rapid  and  highly  sensitive  electrical  measurement  technique  is 
proposed  and  demonstrated  [59]  which  is  designed  to  separate  the  effects  on  the  CV 
and  IV  characteristics  from  stress-induced  areal  nonuniform  interface  and  oxide  traps. 

This  new  technique  combines  the  DCIV  and  AS  methods  to  analyze  samples 
stressed  by  uniform  and  nonuniform  SHEi,  SHEi-f-CHEi,  and  FNTEi  to  uniformly  and 
nonuniformly  charge  and  discharge  the  oxide  traps  via  impact  emission,  low  field 
thermal  capture,  and  EFSE  emission  transitions  discussed  in  chapter  4.  For  example,  it 


89 


90 

was  suggested  that  the  constant  Ig  (unchanged  during  the  low  field  SHEi  stress)  and  the 
changing  distortion  of  the  subthreshold  swing  in  Figure  4.20c  implied  the  presence  of 
areally  inhomogeneous  oxide  charge  generated  during  the  8MV/cm  FNTEi  stress 
(Fig. 4. 20a).  The  design  of  experiments  in  this  chapter  is  based  solely  on  the 
assumption  that  interface  traps  generated  during  the  high  field  (FNTEi  or  SHEi/CHE) 
or  low  field  (SHEi)  phase- 1  will  not  be  removed  or  annihilated  during  the  EFSE 
emission  or  the  low  field  SHEi  neutralization  of  Qqj  during  phase-2.  Therefore  any 
stress-generated  Alg  increase,  due  to  ADjj,  during  phase- 1  will  be  preserved  during 
phase-2.  while  the  stress-generated  areally  nonuniform  AQqj  in  phase- 1  will  be 
neutralized  during  phase-2.  The  results  presented  in  this  chapter  demonstrate:  1)  the 
separation  of  the  generated  interface  traps,  ADfj.  from  the  areally  nonuniform  charging 
and  uniform  neutralization  of  the  oxide  traps,  AQqt-  2)  that  charging  +Qot  by  high 
field  SHEi  with  Vq5=0  is  areally  uniform  whereas  3)  charging  +Qot  by  FNTEi  is 
areally  nonuniform,  and  4)  that  low  field  (<lMV/cm)  -Qqt  charging  at  77°K  is  also 
areally  nonuniform. 

5.1  Nonuniform  +0oj  via  FNTEi 

One  possible  explanation  to  describe  the  nonuniformity  of  the  FNTEi  in 
MOST's  is  that  the  oxide  voltage  drop  (Vq^)  in  the  gate/source  and  gate/drain  overlap 
regions  differs  from  that  in  the  channel  or  well  region  because  of  the  difference  in 
doping  impurity  concentrations,  n+  in  the  source  and  drain,  and  p  or  p-  in  the  channel 
or  well.  It  is  pointed  out  (chapter  3  of  [10,54])  that  the  difference  between  the  gate  bias 
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(Vq)  and  the  oxide  voltage  drop  (Vqx)  is  due  to  three  factors:  1)  the  total  potential  drop 
in  the  polysilicon  gate,  V^.g^^^^,  and  silicon  substrate,  V^.j^^s'  2)  the  gate/substrate  work 
function  difference,  <J>gX'  ^^^  ■^)  ^h^  charged  oxide  and  interface  traps: 

Vox  =  Vg  +  Vs.ga^e  +  Vs-subs  +  ^GX 

=  Vg  +  Vs.g,^,  +  Vs.,„bs  +  VpB  5.1 

where  Oqx  is  equal  to  the  flatband  voltage  (Vpg)  if  Qqt  and  Q,-p  are  negligible  ('X' 
denotes  the  silicon  substrate,  G  the  gate,  and  S  the  gate/oxide  and  oxide/substrate 
interface  position).  For  example,  Vpg  for  a  sample  with  an  n"^-poly-gate  concentration 


of  Ngg  =  1x10^°  phosphorus/cm^  and  a  p-well  concentration  of  Pxx  =  3x10 


16 


boron/cm^  is  Ypg  =  Vp.g,.,  -  Vp.,,,„„„  =  -  (kbT/q)log,(NGGPxxK)  =  -  0.93V  whereas 
VpB  in  the  gate  overlap  regions  of  the  self-aligned  n"^  source/drain  is  Vpg  =  - 
(kBT/q)log,(NGG/NDD)  =  OV.  In  addition  to  the  difference  in  Vpg,  the  silicon  surface 
energy  band  bending  in  the  p-well,  V^.^^,,  or  V^.^hannd'  will  be  greater  than  that  in  the 
n-'source/drain  regions,  Vs.,„„,^^  and  Vs.d.^i^  because  of  the  lower  dopant  impurity 
concentration  in  the  well  and  channel  and  higher  dopant  impurity  concentration  in  the 
drain  and  source,  in  addition  to  their  opposite  conductivity  type.  The  total  potential 
drop  in  the  silicon  gate  and  substrate,  Vg.g^,,,  Vg.^hannei'  Vs-drain  and  Vs.,„u,^e'  has  two 
contributions:  1)  the  potential  drop  in  the  accumulation  or  inversion  layer  (Vg.^^^  or 
Vs.j„y)  and  2)  the  potential  drop  in  the  surface  space  charge  layer  when  biased  into 
inversion.  The  second  component  is  zero  in  accumulation  and  twice  the  bulk  potential 
(2Vp)  in  inversion.  Hence,  given  a  particular  gate  bias  condition,  Vg,  and  neglecting 
any  potential  drop  in  the  highly  doped  polysilicon  gate,  the  oxide  voltage  drop  in  the 
region  over  the  p-weU  and  the  source/drain  overlap  regions  is 
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Vox-wei,  =  Vg  -  VpB  -  Vs.„,n  5.2 

^OX-source/drain  ~  ^G  ~   *  FB  ~  ^S-source/drain  ~  ^G  ^-^ 

Therefore  during  the  accumulation  stress  of  an  n-BiMOST  by  FNTEi  from  the  poly- 
gate  (Vq  =  -13.6V),  the  voltage  drop  in  the  oxide,  Vox»  will  be  higher  at  the  edge  of  the 
channel  and  in  the  source  and  drain  overlap  regions  than  in  the  the  channel  over  the  p- 
well  since  IVpgl»IVs.^.^„l  or  IVqxI  =  IVcl-IVpsl.  The  higher  Vqx  in  the  source/drain 
regions  would  then  give  rise  to  significantly  higher  Fowler-Nordheim  tunneling 
injection  current  passing  through  the  edge  and  source-drain  regions  of  the  oxide  due  to 
the  exponential  dependence  of  Iq  on  Eqx=Vqx/Xo.  Furthermore,  the  higher  Vqx  near 
the  source  and  drain  will  result  in  increased  -i-Qot  generation  in  these  regions  since 
more  carriers  will  gain  the  required  kinetic  energy  (>7eV)  to  impact  emit  the  trapped 
electrons  at  the  oxygen  vacancy  sites.  The  results  of  this  case  were  presented  in 
Figures  4.16-4.20.  Figure  4.20c  (phase-2)  illustrated  the  recovery  in  AS  but  not  %AIg 
during  the  low  field  SHEi  which  neutralized  the  nonuniform  -i-Qot- 

Conversely,  under  the  inversion  stress  (Vq  >  0  in  nMOST)  the  surface  in  the  p- 
well  region  is  inverted,  and  electrons  in  the  silicon  substrate  will  tunnel  through  the 
gate  oxide  layer  into  the  poly-Si  gate.  The  negative  flatband  voltage  over  the  p-well 
would  increase  Vq^  by  approximately  0.93V  in  the  center  of  the  channel  over  the  p- 
well  compared  to  the  source  and  drain  edge  regions.  However,  the  flatband  voltage  is 
now  compensated  by  the  potential  drop  in  the  p-silicon  well  or  substrate  to  invert  the 
surface  to  n-type  and  create  the  n-channel.  This  potential  drop  is  Vs.^.^^^^g,.i„^  =  2Vp.^g,j 
so  the  total  voltage  difference  between  V^  and  Vq^  is  rather  small,  2Vp.j,^3^jjg[  + 
Vc  .      ,  •    >  =0. 

S-channel-mv  —      "• 
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Thus,  a  FNTEi/SHEi  stress-and-measure  or  SAM  experiment  was  conducted 
with  positive  V,^  applied  polarity  to  the  poly-gate  (inversion)  in  order  to  tunnel 
electrons  from  the  silicon  substrate  into  the  oxide  and  to  impact  emit  the  trapped 
electrons  at  the  neutral  nonbridging  oxygen  vacancy  centers  in  order  to  buildup 
positive  oxide  charge.  The  source,  drain,  well  and  emitter  of  the  BiMOST  were  all 
grounded  with  respect  to  the  gate  electrode  which  was  set  to  Vq  =  +13. 6V  to  determine 
if  the  areally  nonuniform  positive  oxide  charge  distribution  could  be  produced  during 
the  inversion  FNTEi  stress  of  this  nBiMOSTs.  Figure  5.1a  illustrates  the  negative  gate 
voltage  shift,  AV^,  during  the  high  field  FNTEi  portion  of  the  stress  and  the 
subsequent  recovery  (Fig.  5.1b)  during  the  low  field  thermal  electron  capture  (SHEi)  to 
neutralize  the  +Qot-  Figure  5.1c  illustrates  the  corresponding  increase  of  the 
subthreshold  swing,  AS,  (measured  at  Io=5nA  and  Vp^^lOOmV)  and  the  percent  shift 
in  the  peak  base  current,  %Al3  (measured  at  Vg£=0.5V)  during  phase- 1.  Figure  5. Id 
illustrates  a  recovery,  by  -60%,  in  AS  towards  its  original  value  while  the  %AIb  due  to 
ADj-p  remains  constant  during  the  phase-2  low  field  SHEi  which  neutralized  the  +Qot- 
This  is  essentially  the  same  result  obtained  under  accumulation  FNTEi  stress  discussed 
in  the  previous  chapter  (Figs.  4.16  and  4.20).  The  recovery  of  AS  indicates  that  an 
areally  nonuniform  distribution  of  positive  oxide  charge  was  generated  during  the  high 
field  FNTEi  which  was  then  neutralized  during  the  low  field  SHEi.  The  coiresponding 
increase  in  Ig  during  the  high  field  FNTEi  without  the  subsequent  recovery  during  the 
low-field  SHEi  is  an  indication  of  the  expected  +Qij  generated  by  high-field  FNTEi 
stress.  These  results  are  consistent  with  the  theory  and  demonstrate  the  separation  of 
Qjj  and  the  areally  nonuniform  +Qq-p. 
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5.2  Room  Temperature  SHEi  Stress 

Positive  charge  can  also  be  generated  by  SHEi  using  the  BiMOST  structure  [56] 
to  give  both  areally  uniform  and  nonuniform  distributions  of  +Qot  ^s  long  as  Vqx 
exceeds  the  7V  threshold.  The  uniformity  of  the  SHEi  can  be  controlled  by  the  d.c. 
voltage  applied  between  the  source  and  the  drain  (Vpj^)  which  alters  the  gate  oxide 
electric  field,  Eqj^,  in  the  y-direction  along  the  channel.  Thus  two  SAM  experiments 
were  performed,  one  with  Vps^OV  during  the  phase- 1  SHEi  to  verify  the  assumed  areal 
uniformity  of  Qqj  and  Dj-p,  and  the  other  with  Vos=4V  during  phase- 1  SHEi  to 
generate  the  nonuniform  Qq-j-  and  Djj.  The  results  provide  further  support  for  the 
methodology  of  separating  Djj  from  areally  nonuniform  -i-Qot- 

5.2.1  Uniform  Oxide  and  hitert'ace  Traps 

The  first  experiment  employs  the  assumed  (and  confirmed  by  the  results) 
uniform  SHEi  stress,  due  to  setting  Vqs=OV.  During  phase- 1,  the  high-field 
(EQx=7.3MV/cm;  Vc=12.8V)  SHEi  stress,  +Qot  ^^^  generated  in  the  oxide  as 
confirmed  by  the  negative  AV^j  (Fig.  5.2a)  via  the  same  impact  emission  process  at 
the  oxygen  vacancy  center  discussed  previously.  This  is  then  followed  by  the  low-field 
(2MV/cm)  SHEi  neutralization  (AVq  recovers)  of  the  +Qot  (Fig.  5.2b).  Both  AS  and 
Alg  increased  during  phase- 1  and  were  roughly  proportional  to  each  other  as  expected. 
Figure  5.2d  shows  that  both  AS  and  Alg  remained  constant  during  phase-2  indicating 
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that  AS  must  all  be  due  to  ADj^  and  not  areally  nonuniform  +QQy.  Note  also  that  Alg 
was  plotted  instead  of  %AIg  in  order  to  illustrate  the  conversion  to  surface 
recombination  velocity  change,  ASq.  The  ratio  AS/AIg  =  (0.50mV/decade)/12.5nA  = 
40kQ/decade  and  ADi/Alg  =  10'°(states/eV-cm2)/12.5nA  =  8xlO*(states/eV-cm2)/nA 
also  illustrates  the  extremely  high  sensitivity  of  the  new  measurement  technique  for 
detecting  interface  traps.  The  change  in  the  surface  recombination  velocity  was 
calculated  using  (4.4)  and  Vgg=0.5V.  Assuming  equal  carrier  capture  cross-sections 
aj=G  =Gq  for  electrons  and  holes  at  the  interface  trap,  the  recombination  cross-section 
is  given  by  (4.5)  which  is  repeated  here 

Cq  =  (2/7t)ASo/(0th  ANjt)  (5.4) 

For  an  interface  trap  energy  distribution  width  of  0.025eV  (about  one  kgT)  and 
0tj,=  lO'^cm/s 

ao  =  (2/7r)ASy(0,h.ADn-.kgT) 
=  (2/7r)(10cm/s)(10^ .  10^°  .  0.025)  =  25k^  =  nr^  =  7t(3A)2  (5.5) 

or  a  radius  of  -SA  which  is  consistent  with  the  bond  length  of  a  localized  dangling 
silicon  bond. 

5.2.2  Nonuniform  Oxide  and  Interface  Traps 

To  complement  the  areally  uniform  +Qot  results  of  the  above  experiment  a 
second  SAM  experiment  was  run  in  which  areal  nonuniformity  in  Qqt  and  Qj^-  was 
generated  during  phase- 1  using  the  SHEi  and  CHE  stress.  This  was  accomplished  by 
applying  a  voltage  between  the  source  and  drain,  Vds=Vd-Vs>0  which  causes  the 
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oxide  voltage  drop  to  vary  from  its  highest  value  at  the  source,  Vq5=Vq-Vqjs-Vs,  to 
its  lowest  value  at  the  drain.  VQp=VQ-Vc-j-p-Vp.  The  ^qts  '^^^  ^gtd  '^^^  ^^e  gate 
threshold  voltages  at  the  source  and  drain  ends  of  the  channel  respectively.  Thus,  if 
V(^-V(3-j-D-Vp  or  the  oxide  voltage  drop  is  less  than  the  threshold  voltage  (~7V)  for 
electron-impact  emission  at  the  Vq  center,  then,  +Qot  could  not  be  generated  in  a 
length  of  the  channel  next  to  the  drain.  Then  +Qq.p=0  in  the  length  of  the  channel  from 
the  drain.  y=yD.  to  y^y;,  where  VQ(yy)=7V.  and  +Q()7  will  increase  from  QQ^(y7)=0  to 
+AQox(y=0)  at  the  source  (y-direction  along  the  length  of  the  channel).  A  sensitive 
and  reliable  d.c.  measurement  of  Qqx  and  Nj^  profiling  (i.e.  their  y-position  or  channel- 
position  dependence)  will  be  reported  [60].  A  highly  peaked  ADj^  or  AQ[y  at  the  drain 
junction  is  not  expected  here  because  Vj)5<Vgs  so  that  the  channel  was  not  depleted 
and  there  was  no  high-field  depleted  space-charge  region  near  the  drain  junction. 
Thus,  the  channel  electrons  were  essentially  at  thermal  equilibrium  in  this  long 
channel.  Hot  electrons  with  high  kinetic  energies  would  be  expected  in  shorter 
channels  when  the  longitudinal  electric  field  over  a  major  part  of  the  channel  length  is 
substantially  larger  than  the  critical  electric  field  for  hot  electrons.  Using  the  linear 
mobility  obtained  in  Fig.4.3a,  the  critical  field  can  be  estimated  from  E^■^^  =  Q^J\i^^^  = 
10'(V/cm)/500(cmVV-s)  =  20kV/cm  =  2V/|im  [10.  pp. 668-669  and  Fig. 3 14.1  on 
p. 252].  Nevertheless,  even  for  long  channels  at  low  longitudinal  electric  fields,  a 
somewhat  nonuniform  distribution  of  AQ^  could  be  expected  from  the  kinetic  energy 
spread  of  hot  holes  which  are  back-injected  from  the  polysilicon-gate  into  the  oxide  as 
recently  demonstrated  by  Lu  and  Sah  [61].  These  oxide  holes  can  drift  to  the 
oxide/silicon  interface  and  break  the  weak  interfacial  hydrogen  and  silicon-oxygen 
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bonds  to  create  a  nonuniform  distribution  of  interface  traps.  The  voltage  conditions 
were  Vg=13.6V,  Vd=4.0V,  Vs=OV  (grounded),  Vg  =  -4V  and  Vbe  =  0.55V.  The 
AVq^,  Alg  and  AS  results  are  shown  in  Figures  5.3a-5.3d.  As  in  the  previous 
experiments  IAVGjl»AS=0.5mV  so  that  AV^x  is  again  mainly  due  to  +AQqj,  but  AS 
should  also  have  a  nonuniformity  component  (as  in  the  FNTEi  case,  Figs.  5.1c  and 
5. Id)  due  to  the  nonuniform  +AQox  and  possibly  nonuniform  AQjj  which  will  now  be 
confirmed. 

The  AVgx  data  (Figs.  5.3a  and  5.3b)  were  obtained  in  strong  inversion  at  Ij)=l|uiA 
(filled  circles  and  triangles)  and  lOpiA  (open  circles  and  triangles).  The  data  labeled 
'Reverse'  in  this  experiment  were  measured  with  the  drain  and  source  interchanged  to 
monitor  the  damage  (i.e.  interface  traps)  near  the  drain  junction  since  the  subthreshold 
Iq-Vq  characteristics  are  sensitive  to  the  interface  traps  located  near  the  injecting 
junction  (the  source  when  Vq=0.1V  or  the  drain  when  Vs=0.1V).  During  phase- 1  the 
near  coincidence  of  the  forward  and  reverse  AVgj  in  strong  inversion  (lO^iA)  and  the 
30%  difference,  AVGT.reverse/^VGT.forward=(-43mV)/(-33mV)=1.3,  at  a  weaker  inversion 
(IjiA)  indicate  that  the  nonuniform  AQqt  and  AQj^-  are  screened  out  by  the  high 
electron  density  in  the  inversion  channel.  This  is  consistent  with  the  results  of  Figs. 
4.19a  and  4.25  which  show  that  the  mobility  is  less  influenced  by  Qox  at  higher  values 
of  inversion  charge  density  Njnv  The  +Qot  charge  buildup  (-AVg^  shift)  during 
phase- 1,  (Fig.  5.3a)  is  completely  neutralized  during  phase-2  as  Yq-^  recovers  (Fig. 
5.3b). 

The  AIb  and  AS  measured  in  the  subthreshold  range  are  shown  in  figures  5.3c 
(phase- 1)  and  5.3d  (phase-2).  The  results  are  similar  to  those  of  the  uniform  stress  in 
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figures  5.2c  and  5. 2d  except  that  the  nonuniform  AQqt  gives  a  large  source-drain 
difference.  AS,„^,^^,  =  ^ASj^^i^.  Figure  5.3d  shows  that  AS,„„^^,,  recovers  during  phase-2 
from  its  maximum  after  phase-1  by  =  0.63  -  0.23  =  0.4mV/decade  due  to 
neutralization  of  the  nonuniform  +Qot  while  the  remaining  ASj-j-  =  0.23mV/decade 
from  ADjT  does  not.  This  is  confirmed  by  the  fact  that  the  ratio 
ADi-j-/AlB=5xlO^(states/eV-cm^)/6.25nA=8xlO^(states/eV-cm^)/nA  remained  the  same 
as  in  the  uniform  SHEi  stress  case  calculated  from  figure  5. 2d. 

5.3    Low  Temperature  SHEi  Results 

In  the  last  chapter  it  was  also  shown  that  significant  areal  nonuniformity  was 
obtained  when  using  low-field  SHEi  (IMV/cm)  to  charge  the  shallow  E^-leV  electron 
traps  in  the  oxide  at  77°K.  Since  high  field  SHEi  at  room  temperature  was  shown  to 
give  uniform  positive  charge  generation  in  section  5.2  the  low-temperature  low-field 
SHEi  experiment  was  repeated  while  measuring  the  top-emitter  base  current  [12]  to 
monitor  any  increase  in  ADjj.  The  top-emitter  Ig  measurement  makes  use  of  the  source 
or  the  drain  as  the  emitter  for  minority  carrier  injection  into  the  base  region  while 
monitoring  the  V^g-Ig  characteristic.  For  this  experiment  the  drain  junction  was  used. 
The  dominant  sources  of  base  current  in  the  MOST  channel  and  base  region  (Fig.  5.4) 
represent  the  following  mechanism: 

^B-i-buik  ""  space  charge  region  recombination  in  the  bulk;  n=2 
iB-iB-buik '"  quasi-neutral  Base  region  recombination  in  the  bulk;  n=l 
^B-iE-buik  "■  quasi-neutral  Emitter  region  recombination  in  the  bulk;  n=l 
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lB-2-surface  "  spacc  choTge  region  recombination  at  the  surface:  n=2 
iB-iB-surface ""  Quasi-neutral  Base  region  recombination  at  the  surface;  n=l 
The  n  represents  the  exponential  dependence  of  the  base  current  (I^  °=  exp[qV/nkgT]) 
where  n=l  is  the  ideal  Shockley  diode  current.  Jj,  and  n=2  is  the  Sah-Noyce-Shockley 
(generation-recombination)  diode  current,  J2  [10].  The  first  three  sources  of  base 
recombination  current  are  from  the  bulk  region  of  the  base  layer  and  hence  not 
modulated  by  V^g  and  therefore  contribute  only  to  the  baseline  Ig  level.  On  the  other 
hand  Ig_2s  and  Ig.igj  will  vary  with  Vcg  to  give  a  peak  in  Ig  at  the  maximum  interface 
recombination  condition.  One  very  important  difference  between  the  bottom-emitter 
and  top-emitter  measurement  is  the  top-emitter's  additional  sensitivity  to  electron-hole 
recombination  at  the  interface  traps  in  the  drain/base  n-i-/p  junction  space  charge  layer, 
Ig  2s,  when  it  is  forward  biased.  In  contrast,  during  the  bottom  emitter  Ig  measurement, 
the  n^/p  drain/well  and  source/well  junctions  are  at  zero  or  reverse  bias  so  that 
recombination  is  negligible  in  the  drain/well  and  source/well  junctions.  Clearly,  the 
gate  bias  at  which  maximum  interface  recombination  occurs  will  differ  between  the 
space  charge  layer  near  the  drain  (Ig_2s)  compared  to  the  quasi-neutral  p-well  (Ib-ibs) 
due  to  the  difference  in  the  surface  potential  of  these  two  regions.  At  low  forward 
junction  bias,  the  Jj  component  will  dominate  since  most  of  the  injected  minority 
carriers  (electrons)  by  the  n-i-drain/p-well  junction  will  recombine  in  the  space-charge 
layer.  Another  key  point  is  that  the  top  emitter  base  current  peak,  due  to  the  lg_2s 
component,  is  highly  sensitive  to  the  area  of  the  space  charge  layer  underneath  the  gate 
oxide  and  will  therefore  be  influenced  by  the  local  doping  impurity  concentration  in 
the  p-well,  such  as  that  of  the  lowly-doped-drain. 
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The  results  indicate  that  as  the  oxide  becomes  negatively  charged  during  the  low 
field  (IMV/cm)  SHEi  stress  (phase- 1)  the  Ib-Vgb  characteristics  shift  in  the  +Wq 
direction  and  le-peak  ^^^^  decreases  significantly  by  almost  40%  (Fig. 5. 5).  The 
decreasing  shift  in  peak  base  current,  Alg,  and  the  corresponding  increase  in  the 
subthreshold  Iq-Vq  swing,  AS,  are  shown  in  Fig.  5.6.  The  +Vq  shift  in  the  Ci^y-V^ 
characteristics  and  the  gate  voltage  shift  for  constant  values  of  Cj^v  -^Vqc,  are  also 
obtained  and  shown  in  Figs.  5.7  and  5.8  respectively,  to  illustrate  the  large  amount  of 
distortion  in  the  CV  during  the  phase- 1  charging.  Next,  the  electrons  trapped  at  the 
shallow  E^-leV  level  are  field-emitted  just  as  in  the  experiment  presented  in  chapter  4. 
The  second  portion  (EFSEe  or  phase-2)  of  Figs.  5.5-5.8  illustrate  a  complete  recovery 
of  the  base  current,  the  subthreshold  swing,  and  the  inversion  layer  capacitance  to  their 
original  prestress  values.  Under  the  assumption  that  interface  traps  generated  during 
phase- 1  would  not  be  annihilated  during  phase-2,  it  would  appear  that  the  distortion  in 
Ig,  AS,  and  Cj^v  are  entirely  due  to  the  areally  nonuniform  negative  charging  of  the 
gate  oxide,  -QotCy)-  The  decrease  in  Ig  during  phase- 1  is  attributed  to  the  nonuniform 
negative  charging  of  the  oxide  near  the  drain/base  junction  which  effectively  increases 
the  concentration  of  holes  near  the  Si02/Si  interface  and  therefore  reduces  the  space 
charge  layer  area  underneath  the  gate  where  electrons  and  holes  may  recombine.  This 
effect  is  illustrated  by  the  region  enclosed  by  the  two  dashed  curves  and  dashed  squares 
in  the  cross-sectional  view  of  the  drain/base  junction  in  Fig.  5.9.  The  nonuniformity 
encountered  during  the  -Qqt  charging  phase- 1  may  be  attributed  to  a  voltage  drop  from 
the  center  of  the  long  channel  (150|im)  towards  the  source  and  drain  (collector) 
junctions  during  the  low  field  SHEi  and  would  be  proportional  to  Rch*Ic-  Such  a 
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Fig.    5.5     The  shift  in  the  base  current  characteristics  during  low  field  (IMV/cm) 
SHEi  stress  to  fill  the  shallow  Ec-leV  electron  traps  (solid  arrows), 
followed  by  EFSE  emission  to  detrap  the  electrons  (dotted  arrows). 
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Fig.    5.6     The  percent  shift  in  the  peak  base  current  and  subthreshold  current  swing 
during  low  field  (IMV/cm)  SHEi  stress  to  fill  the  shallow  Ec-leV  electron 
traps  (left  of  dashed  line)  followed  by  EFSE  emission  to  detrap  the  electrons 
(right  side  of  dashed  line). 
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Fig.    5.7     The  shift  in  the  lOkHz  inversion  capacitance,  Cj^^  vs  d.c.  gate  voltage,  Vg, 
during  low  field  (IMV/cm)  SHEi  stress  to  fill  the  shallow  Ec-leV  electron 
traps  ^direction  of  solid  lines)  followed  by  EFSE  emission  to  detrap  the 
electrons  (direction  of  dashed  Hnes). 
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Fig.    5.8     The  distortion  in  the  lOkHz  inversion  capacitance,  C^^^  vs  d.c.  gate  voltage 
during  low  field  (IMV/cm)  SHEi  stress  to  fill  the  shallow  Ec-leV  electron 
traps  (direction  of  sohd  arrow)  followed  by  EFSE  emission  to  detrap  the 
electrons  (direction  of  dashed  dashed). 
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voltage  drop  would  significantly  alter  the  accelerating  electric  field  profile  in  the 
silicon,  E^,  and  hence  result  in  an  areally  nonuniform  injection  profile. 

In  this  chapter  a  new  and  highly  sensitive  methodology  for  detecting  and 
separating  areal  nonuniformity  in  the  oxide  charge  distribution  from  the  stress- 
generated  interface  traps  was  demonstrated  [59].  The  methodology  employed  the  two 
measurements:  1)  the  new  DCIV  method  recently  proposed  by  Neugroschel  and  Sah 
[12]  to  detect  the  stress-generated  interface  traps  by  measuring  the  Ig  of  the  BiMOST, 
and  2)  the  traditional  subthreshold  dc  Iq-Vq  swing  to  detect  both  the  stress-generated 
interface  traps  and  nonuniform  oxide  charge  distributions.  The  methodology  was 
demonstrated  in  five  SAM  experiments  performed  on  five  BiMOSTs  where  spatially 
uniform  and  nonuniform  (from  source  to  drain)  oxide  charge  distributions  were 
generated  by  experiments  using  the  two  FNTEi  polarities,  in  two  SHEi/CHEi 
experiments  with  Vps^O  and  Y^^=4Y,  and  by  low  field  SHEi  at  77°K.  The  increased 
distortion  of  the  a.c.  C^^^-Yq  characteristics  during  the  low  field  SHEi  for  -Qqt 
charging  and  its  complete  recovery  provided  further  support  to  the  new  combined  d.c. 
AS  and  DCIV  measurement  methodology  proposed  in  this  chapter. 


CHAPTER  6 

CURRENT-ACCELERATED  CHE  STRESS 

FOR  RAPID  TIME-TO-FAILURE  EXTRAPOLATION 


Thus  far  the  injection  stress  applied  to  the  n-BiMOST's  in  this  study  included 
FNTEi,  SHEi  and  SHEi/CHE.  This  chapter  compares  the  results  of  the  conventional 
voltage-accelerated  CHE  and  the  new  current-accelerated  CHE  stressing 
methodologies  and  their  effect  on  the  channel  mobility,  the  drain  saturation  current,  the 
gate  threshold  voltage  and  the  DCIV  peak  base  current.  The  traditional  time-to-failure 
(TTF)  extrapolation  by  voltage-accelerated  CHC  stress  will  also  be  reviewed  to 
illustrate  its  limitations  and  to  point  out  the  advantages  of  the  current-accelerated 
methodology. 

6.1  Conventional  Channel  Hot  Carrier  Stress 

The  voltage-accelerated  channel  hot  carrier  (CHC)  stress  methodology,  for  n- 
channel  and  p-channel  MOST's  has  become  an  industry  standard  for  extracting  the 
time-to-failure  (TTF)  of  MOST's  to  assess  the  reliability  of  new  technologies  at  lower 
operation  voltages.  The  customary  CHC  stress  condition  is  to  apply  Wj^=2.5Vq,  so  that 
the  channel  is  strongly  depleted  at  the  drain,  and  the  substrate  (the  well  in  MOST's  or 
base  in  the  BiMOST)  current  is  maximized  during  stress.  To  avoid  confusion  with  the 
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base  current  measurement  (Ib-Vgb)  used  to  characterize  interface  states  in  the  previous 
two  chapters  the  substrate  current  will  be  denoted  by  the  standard  notation  I^  (or  I^yg) 
although  they  are  the  very  same  current  measured  at  the  same  terminal,  B,  X,  or  W 
(Well),  but  at  different  bias  configurations.  The  substrate  current  component  originates 
from  the  impact  generation  of  electron-hole  pairs  by  energetic  electrons  transiting 
through  the  high-field  space-charge  region  of  the  reverse  biased  drain/base  junction 
(Fig.  6.1).  The  inversion  electrons  entering  the  channel  at  the  source  drift  through  the 
inversion  channel  until  they  reach  the  drain-junction  space-charge  layer  where  they  are 
then  accelerated  by  the  high  electric  field  to  high  kinetic  energies.  If  the  electron 
kinetic  energy  gained  is  more  than  1.2eV,  as  recently  demonstrated  theoretically  and 
experimentally  by  Lu  and  Sah  [62],  before  exiting  the  drain,  the  electrons  can  impact 
release  an  electron-hole  pair.  The  negative  sign  of  I^  arises  from  the  impact-generated 
holes  exiting  through  the  p-well  contact.  Furthermore,  if  these  energetic  electrons  gain 
enough  kinetic  energy  (3.1eV)  they  may  also  be  injected  into  the  oxide  to  give  negative 
[56,57]  or  positive  [56,  61,  62]  charge  as  demonstrated  in  the  preceding  chapters,  or  to 
break  the  weak  interfacial  bonds,  such  as  the  Si-H  and  SiO-H  discussed  by  Sah  [9],  and 
create  an  interface  state,  Qj-p. 

Typical  Ix-Vqs  characteristics,  with  V^s  as  the  parameter,  taken  on  the 
100xlOO|im^  n-channel  BiMOSTs  used  in  this  study,  are  shown  in  Fig.  6.2.  These 
conventional  CHE  curves  illustrate  the  limiting  nature  of  this  technique  when 
extrapolating  TTF  at  low  operating  voltages  {W^  and  Vg)  within  a  reasonable 
measurement  time  of  150  hours  or  less  (<5  days).  The  conventional  approach  therefore 
is  to  accelerate  the  stress  by  raising  the  voltages  and  corresponding  kinetic  energy  of 
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Figure  6.2      Substrate  hole  current,  Ix,  versus  gate  bias  at  Vd=5. 5-6.5. 
The  dotted  lines  represent  bottom-emitter  current- 
accelerated  data  with  Vbe=0.55V 
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the  carriers  responsible  for  the  damage  at  the  Si02/Si  interface.  In  doing  so  it  is 
incorrectly  assumed  that  the  microscopic  failure  mechanisms  are  independent  of  the 
spatial  and  kinetic  energy  distributions  of  the  carriers  (electrons  in  the  nMOST). 
Clearly  therefore  the  traditional  voltage-acceleration  method  is  not  reliable  for 
extrapolating  TTF  in  future  deep  submicron  technologies  operating  at  voltages  below 
2V  as  will  be  shown  by  the  results  in  the  next  section. 

6.2  Bottom-Emitter  Current- Accelerated  CHE 

Recently  Neugroschel  and  Sah  proposed  two  current  acceleration  techniques 
[12,13]  for  extrapolating  the  operation  TTF  at  low  voltages  in  BIT  with  stress-time 
acceleration  greater  than  100.  In  this  thesis  a  similar  current- acceleration  method  for 
CHE  TTF  extrapolation  is  demonstrated  for  the  first  time  on  nMOSTs.  The  current- 
acceleration  for  MOSTs  exploits  the  unique  vertical  bipolar  structure  of  the  BiMOST 
test  structure  by  forward  biasing  the  bottom-emitter/base  junction.  The  forward  biased 
minority  carriers  (electrons  in  nMOST  of  this  study)  diffuse  through  the  base  layer 
towards  the  SiOj/Si  surface  where  they  are  swept  across  the  channel  inversion  layer 
along  with  the  carriers  originating  from  the  source  junction  (Fig.  6.3).  Hence  the  stress 
due  to  the  impact  generation  of  electron-hole  pairs  occurring  near  the  drain  is  current- 
accelerated  by  the  larger  electron  flux  that  enters  the  high-field  space-chai'ge  layer  near 
the  Si02/Si  surface.  More  important,  the  kinetic  energy  distribution  of  these  additional 
electrons  is  the  same  as  those  in  the  nonaccelerated  CHC  case  since  the  V^,  and  Vq  are 
kept  constant.  The  alteration  of  the  channel  accelerating  electric  field  distribution  by 
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the  injected  electrons  is  expected  to  be  small.  It  is  noted  however  that  because  of  the 
diffused  emitter  geometry  of  this  n-BiMOST  structure  there  is  a  significant 
contribution  to  the  drain  current  from  electrons  entering  the  drain  through  its  bottom- 
facing  junction  and  entering  laterally  near  the  surface  outside  the  channel. 
Furthermore,  there  will  be  additional  impact  generation  of  electron-hole  pairs 
occurring  in  the  bulk  (bottom-facing  part)  of  the  drain/base  junction  space-charge  layer 
and  near  the  surface  outside  the  channel.  These  secondary  sources  of  drain  current  will 
increase  Ix  during  stress  but  will  not  damage  (or  create  interface  traps  at)  the  Si02/Si 
interface  in  the  channel.  The  secondary  surface  impact  generation  current  component 
outside  the  channel  will  however  increase  the  baseline  of  the  Ig  current  (DCIV)  by 
creating  interface  traps  in  the  oxide  layer  adjacent  to  the  drain.  However,  the  surface 
potential  in  the  location  of  these  interface  states  will  not  be  modulated  by  the  gate  bias 
when  measuring  the  Ib-Vgb  characteristic  and  therefore  will  not  exhibit  the  peak 
behavior  like  the  interface  states  in  the  channel.  The  increase  in  I^  when  Vgg^O  is 
easily  demonstrated  by  the  Ids'^ds  family  of  curves  with  y^^^=\-3V  and  VBg=0-0.55V 
as  secondary  variables  in  Fig.  6.4.  The  value  of  the  accelerating  emitter  current  is  also 
shown  in  Fig.6.4. 

The  change  in  the  top-emitter  (drain)  DCIV  characteristics,  Alg,  during  CHE 
stress  (Vd=8V,  Vg=2.8V  and  VgE^O.OV)  for  the  lOOxlOOiim^  n-BiMOST  are  shown  in 
Fig.  6.5.  The  value  of  the  accelerating  reverse  bias  voltage  on  the  drain/base  junction 
is  Vds  -  Vds,,^  =  Vds  -  (Vgs  -  Vgt)  =  8-(2.8-0.5)  =  5.7V.  Though  the  peak  near  the 
flatband  voltage  or  V^g  =  -  0.38V  is  unaffected  during  the  CHE  stress,  an  increasing 
secondary  Ig  peak  is  observed  in  the  accumulation  portion  of  the  Ig- V^g  curve. 
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Figure  6.4     Drain  current  vs  drain  voltage  for  Vg=1,  2,  3  volts.  The 
solid  marker  curves  represent  the  current-accelerated  data 
with  Ig  (solid  line  without  marker)  plotted  on  the  right  axis. 
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Figure  6.5      The  top-emitter  DCIV  characteristics  versus  CHE  stress  for  an 
n-BiMOST.    The  stress  conditions  were  Vd=8V  and  Vg  =2.8V. 
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indicating  a  large  build-up  of  interface  states  near  the  drain  junction.  The  percent 
increase  in  Ig  of  this  secondary  peak  is  plotted  in  Fig.  6.6  versus  the  drain  fluence,  Nji^j 
=  (l/qZ)jl£,dt.  The  corresponding  gate  voltage  shift  of  the  secondary  Ig  peak  is  plotted 
in  Fig.  6.7  illustrating  a  large  -0.7V  shift  in  the  local  threshold  voltage  near  the 
drain/base  junction  due  to  the  build-up  of  Qu  and  possibly  +Qot-  It  is  also  evident  that 
the  microscopic  failure  mechanisms  responsible  for  the  %AIg  are  limited  as  indicated 
by  the  saturating  nature  of  both  the  %AIb  and  the  AVQ.jgpj.^^  curves.  Hence  the 
empirical  and  statistical  linear  or  log-log  TTF  extrapolation  methodology  traditionally 
used  in  industry  for  estimating  the  damage  at  long  CHE  stress  times  is  clearly  in  error. 

Next  the  drain  bias  was  kept  constant  (Vd=8.0V)  while  the  gate  bias  was  varied 
between  Vq  =  2.8V  and  4.1V  so  that  the  accelerating  bias  near  the  drain/base  junction 
was  varied  between  Vps-V^s-sat  =  5.7V  to  4.4V.  The  increase  in  %AIg  vs  Nji^j  for  the 
three  V^  values  with  V^  =  8.0V  is  plotted  in  Fig.  6.8  showing  a  large  variation  in  the 
saturation  values.  As  the  reverse  bias  field  near  the  drain  decreases  so  does  the  damage 
or  interface  trap  generation  as  exhibited  by  the  %AIb  data  which  would  indicate  that 
there  is  a  threshold  voltage  beyond  which  the  damage  becomes  negligible.  Hence  the 
empirical-statistical  linear  or  log-log  TTF  extrapolation  to  lower  operating  voltages 
traditionally  used  in  industry  is  also  in  error. 

A  second  set  of  CHE  stress  and  measure  (SAM)  experiments  were  conducted  for 
low  Vq  =  0.9V  and  Vp  =  5V  to  8V,  with  and  without  the  current-acceleration  from  the 
forward-biased  bottom-emitter.  The  Alg  vs  stress  time  for  V^,  =  0.9V,  V^  =  6.0V  and 
varying  degrees  of  current-acceleration,  as  controlled  by  the  current-limited  forward 
bias  applied  to  the  bottom  emitter/base  junction,  Vg^  =  O.OV  to  6.0V,  are  shown  in  Fig. 
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Figure  6.6      Percent  increase  in  the  top-emitter  DCFV  secondaiy 

iB-peak  value  versus  the  drain  fluence  during  CHE  stress. 
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Figure  6.7      Gate  voltage  shift  in  the  peak  value  of  the  base  current 
measured  by  the  top-emitter  DCIV  method. 
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Figure  6.8      Percent  increase  of  the  top-emitter  DCIV  secondary 
Is-peak  versus  the  drain  fluence  during  CHE  stress. 
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6.9a.  As  indicated  in  this  figure,  for  increasing  current-acceleration  (increasing  Ygg  or 
Ig)  the  Alg  degradation  is  accelerated  by  as  much  as  100  times.  When  increasing  Vgg 
between  2.5V  and  6.0V  it  appears  that  the  degradation  may  have  actually  decreased 
slightly.  This  can  be  attributed  to  an  increase  in  the  baseline  of  the  Ig  current  in  the 
surface  outside  the  channel  region  as  discussed  earlier  in  this  section.  If  the  data  is 
adjusted  to  account  for  this  increase  in  the  baseline  Ig  the  two  curves  coincide  quite 
well  as  indicated  in  Fig.  6.9b.  The  remainder  of  the  data  presented  in  this  chapter  will 
therefore  be  adjusted  to  account  for  this  effect. 

Nevertheless,  it  appears  that  the  current-acceleration  due  to  the  additional  Iq 
contribution  from  the  bottom-emitter  is  limited  beyond  VgE=2.5  or  lp=2.9mA  possibly 
due  to  current  crowding  at  the  surface.  Any  additional  current  flowing  through  the 
drain  beyond  this  value  (2.9mA)  must  therefore  be  collected  by  the  peripheral  bulk 
area  of  the  drain/base  junction  as  opposed  to  the  inversion  channel  at  the  Si02/Si 
interface.  As  mentioned  earlier  this  current  will  not  contribute  to  the  Alg  damage 
except  by  increasing  the  baseline  Ig  since  it  is  not  flowing  along  the  channel  area.  This 
is  evident  in  Fig.  6.10  which  replots  the  Alg  data  of  Fig.  6.9b  versus  the  drain  tluence. 
The  non-accelerated  data  has  now  merged  with  the  current-accelerated  data  for  VgE=0 
to  2.5V  which  supports  the  theory  that  the  kinetic  energy  of  the  carriers  in  the  cuiTent- 
accelerated  case  has  not  changed  as  would  have  been  the  case  using  voltage- 
acceleration.  The  high  VBg=6.0V  case  does  not  merge  with  the  other  three  lower 
forward-bias  curves  because  the  drain  fluence  has  the  additional  contribution  from  the 
peripheral  bulk  drain  current  just  mentioned.  The  merging  of  the  curves  in  Fig.  6. 10  is 
also  consistent  with  the  fact  that  the  drain  current  as  well  as  the  substrate  current  ratios 
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Increase  of  the  top-emitter  DCIV  secondary  iB.peak 
versus  the  total  injection  time  during  CHE  stress.with 
varying  degrees  of  bottom-emitter  current-acceleration. 
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Figure  6.9b      Increase  of  the  top-emitter  DCIV  secondary  Is-peak 

versus  the  total  injection  time  during  CHE  stress. with 
varying  degrees  of  bottom-emitter  current-acceleration. 
The  data  has  been  adjusted  by  the  baseline  Ig  value.. 
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Figure  6. 10      Increase  of  the  top-emitter  DCIV  secondary  Is-peak 

versus  the  total  drain  fluence  during  CHE  stress. with 
varying  degrees  of  bottom-emitter  cunent-acceleration. 
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before  and  after  acceleration  are  proportional  to  the  amount  of  the  acceleration.  For 
example  the  acceleration  factor  between  Vgg  =  OV  to  0.9V  is  approximately  20  while 
h-^Jh  =  280^A/14^iA  =  20  and  Ix-acA  =  14^A/0.75nA  =  19;  where  Ix.^cc  and  Ix  are 
the  accelerated  and  non-accelerated  values.  The  same  good  agreement  is  also  found  for 
the  8-10  times  acceleration  range  between  VgE  =  0.9  to  2.5V  where  lD.ac</lD  ~  ^^-^  and 
Ix-acc/^x  -  ^■^-  However,  the  excess  accelerating  current  flow  will  modulate  the  base 
carrier  concentration  and  produce  a  2-D  resistive  voltage  drop  in  the  base  region  which 
in  turn  will  alter  the  field  profile  near  the  drain/base  and  the  source/base  junctions. 
Collectively  these  effects  are  evidenced  by  the  -AV^j  of  the  device  (Fig.  6.4)  which 
may  also  account  for  the  divergent  or  non-parallel  shift  of  the  Vg£=2.5V  and  VgE=6.0V 
stress  curves  at  high  tjj^j  (Fig.  6.9b)  since  a  -AV^-j-  would  reduce  the  value  of  V^.^^p  = 
Vps  -  (Vgs  -  Vqj)  and  correspondingly  the  field  near  the  drain/base  junction. 

Now  that  the  functionality  of  the  bottom  emitter  current-acceleration  method  in 
MOST's  has  been  illustrated  the  TTF  and  QTF  (charge  to  failure)  will  be  extracted  for 
Vq  =  0.9V,  VgE  =  0  and  2.5V  and  a  range  of  V^  =  5V  to  8V.  The  charge  to  failure  is 
calculated  from  the  drain  fluence  measured  in  these  experiments.  Since  the  AYqj, 
Al^.^aj,  A^i„^ff  degradation  for  these  large  lOOxlOO^im^  devices  was  too  small  (<0.1% 
change)  the  100%  Alg  value  was  used  as  the  TTF  and  QTF  criteria  in  Figs.  6.1 1  and 
6.12  respectively.  The  TTF  data  illustrate  a  varying  amount  of  acceleration  which 
appears  to  increase  at  lower  values  of  V^.  The  QTF  data  on  the  other  hand  merges  at 
the  lower  V^  stress  voltages  but  deviates  for  V^  >  6V  or  Vu-^u-s^i  =  5.6V.  One  reason 
for  the  deviation  at  higher  voltages  may  be  due  to  the  introduction  of  another  failure 
mechanism  because  of  the  two  dimensional  nature  of  the  injection  near  the  drain  when 
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Figure  6. 11      Charge  to  failure  measured  at  100%AIb  during  CHE 

stress  with  Vg=O.OV  and  Vbe=0  for  the  non-accelerated 
data  and  Vbe=2.5V  for  the  accelerated  data. 
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Figure  6.12     Time  to  failure  measured  at  100%AIb  during  CHE 

stress  with  Vg=O.OV  and  Vbe=0  for  the  non-accelerated 
data  and  Vbe=2-5V  for  the  accelerated  data. 
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using  the  bottom  emitter  to  accelerate  the  stress.  In  fact  at  high  enough  reverse  bias 
across  the  drain/base  junction  it  would  not  be  unlikely  to  have  SHEi  occurring  neai*  the 
drain.  Another  factor  which  may  be  responsible  is  the  spatial  distribution  of  Q[j  or 
Qox  near  the  drain  since  the  top-emitter  DCIV  measurement  was  shown  to  be  sensitive 
to  areal  inhomogeneity  in  addition  to  ADj^  in  chapter  5. 

6.3  Top-Emitter  Current- Accelerated  CHE 

An  alternative  approach  to  the  bottom-emitter  current-acceleration  method  just 
demonstrated  for  long-channel  nBiMOSTs  is  to  forward  bias  the  source  junction 
during  the  CHE  stress,  as  suggested  by  Sah.  The  results  of  this  measurement,  shown  in 
Fig.  6.13,  exhibit  a  very  limited  amount  of  acceleration  when  Vg^  =  -l.OV,  Vq^  = 
0.9V  and  V^is  =  ^^  compared  to  the  previous  condition  of  V^x  =  OV,  V^x  =  0.9V  and 
Vj5s  -  7V.  The  limited  acceleration  in  this  case  is  most  likely  due  to  the  lower  voltage 
drop  in  the  depleted  drain  space-charge  layer  since  V^^  =  V^x  -  V^x  =  0.9V  -  (-1  .OV) 
=  1.9V  so  that 

Vo-depCaccelerated)  =  V^s  -  (V^s  -  V^t) 
=  7V-(1.9V-Vgt.o  +  AVgt) 
=  5.1V  +  Vgt.o-AVgt 
while  the  non-accelerated  case  has 

VD-dep(non-accelerated)  =  V^s  -  (V^s  -  VpT-o) 
=  7V-(0.9V-Vgt.o) 
=  6.1V  +  Vgt.o- 
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Figure  6. 13      Shift  in  the  base  current,  Ib,  due  to  nonaccelerated  and 
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top-emitter  (source)  current  accelerated  CHE  stress  at  equal  Vp. 


137 

Indeed  by  comparing  this  source  current-accelerated  CHE  stress  case  to  the  SAM 
results  at  Vsx=OV,  Vgx=0-9V  and  V£,s=6V  a  large  acceleration  is  observed  as  indicated 
in  Fig. 6. 14.  In  addition,  proximity  of  the  bottom-emitter  junction,  which  acts  as  the 
collector  in  this  forward-biased  source  configuration,  may  actually  be  forward  biased 
near  the  drain  if  there  is  a  significant  lateral  (y-direction)  current  flowing  in  the  base  to 
cause  a  lateral  base  voltage  drop.  In  any  case,  this  current-acceleration  scheme  is  most 
desirable  for  the  short-channel  devices  since  the  limited  collection  area  of  the  short 
channel  would  limit  the  bottom-emitter  acceleration  efficiency  and  the  encroachment 
of  the  drain  and  source  space-charge  regions  would  reduce  the  bottom-emitter  injected 
electrons  that  reach  the  channel  as  indicated  by  the  cross-sectional  view  in  Fig.6.3. 

A  final  set  of  CHE  stress  measurements  were  conducted  on  short-channel 
(L=1.6|im,  W=100|im)  nBiMOSTs  to  illustrate  the  current-acceleration  on  AV^.j^t, 
AIq.j^,,  An„,jjj,  and  AV^.^^  since  these  parameters  are  directly  correlated  to  the  device 
performance.  In  addition,  a  first  order  correction  in  V^x  was  made  to  compensate  for 
the  shift  in  the  threshold  voltage,  -AVqj,  due  to  the  forward  bias  applied  on  the 
source/well  junction,  Vsx=-0.5V.  The  results  of  the  CHE  stress  in  the  aforementioned 
parameters  are  shown  in  Figs.  6.15-6.18.  The  nonaccelerated  stress  was  conducted  at 
Vsx=OV,  Vgx=6.0V  and  Vdx=2.0V  (Vex=0)  while  the  source-emitter  current 
accelerated  stress  had  Vsx=-0.5V,  Vgx=2.0V  and  Vdx=6.2V.  It  is  noted  that  the 
degradation  in  the  drain  saturation  current,  AI^..^^  extracted  at  constant  values  of 
Vd=Vg  (Fig.  6-16)  is  primarily  due  to  the  degradation  in  the  effective  saturation 
mobility,  Aii^^^^j  (Fig.  6-17).  An  acceleration  factor  of  ~8  is  illustrated  for  AIq.,^,,  A\x„^^^, 
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Figure    6.17      Shift  in  the  saturation  effective  mobility  from  a  fit  to  the  Sah- 
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reverse  I-V  data  the  source  and  drain  have  been  interchanged. 
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and  AVg_,j„,  extracted  at  Ip=100^A,  while  in  contrast  the  corresponding  gate  threshold 
voltage  shift,  AVgj_j,^j,  extracted  in  the  saturation  range  was  <lmV.  The  discrepancy 
between  AVci-.^at  and  AV^.^^  is  due  to  the  fact  that  the  dominant  Afi^j^j  contribution  was 
factored  out  of  the  degradation  by  the  two  parameter  (VcT-sat-  M^nsat)  fit  of  the  saturation 
data  to  the  Sah-Pao  theory  [55]. 


CHAPTER  7 
CONCLUSIONS 

Novel  measurement  techniques  and  results  are  obtained  in  this  thesis  and 
described  in  chapters  2-6.  The  primary  focus  of  this  thesis  work  is  to  investigate  and 
model  the  degradation  mechanisms  of  the  current  and  mobility  of  the  surface- 
inversion-channel  silicon  MOS  transistor  for  improving  the  design  methodology  of 
future  generations  of  deep  sub-half-micron  integrated  circuits.  Current  and  mobility 
degradations  in  the  linear  and  current-saturation  ranges  of  the  MOS  transistor  current- 
voltage  characteristics  are  measured  under  controlled  and  reversible  oxide  charging 
and  discharging  stress  conditions  to  simulate  the  electrical  stress  experienced  by  MOS 
transistors  during  operation.  The  controlled  stresses  employed  in  the  stress-and- 
measure  (SAM)  experiments  of  this  investigation  included:  (1)  substrate  hot  electron 
injection,  (2)  Fowler-Nordheim  electron  tunneling  injection,  and  (3)  channel  hot 
electron  stress.  The  measurements  included:  (1)  the  linear  drain  current  versus  gate 
voltage,  Id-Vq,  (2)  the  saturation  lo-sat-Vc-sat'  (3)  the  inversion  layer  capacitance  versus 
gate  voltage,  C-^^-Ygx>  and  (4)  the  base  current  versus  gate  voltage,  Ib-Vgb'  with 
forward-biased  top-emitter  or  bottom-emitter  base  current,  via  the  new  DCIV  method 
recently  invented  by  Neugroschel  and  Sah  and  demonstrated  by  them  and  their 
graduate  students. 

An  in-depth  literature  review  is  given  in  chapter  2.  The  various  physical 
mechanisms  affecting  the  carrier  (electron  or  hole)  mobility  in  the  semiconductor 
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surface  channel  investigated  by  previous  authors  are  described  and  analyzed  or 
summarized. 

Chapter  4  describes  the  effects  of  varying  the  oxide  electric  field  (4MV/cm  and 
5MV/cm)  during  SHEi  stress  on  the  degradation  rate  of  AI^.^^,,  AVgt^  A|lI,,„.  A|I,^,.  AS 
and  Alg.  Two  groups  of  novel  experiments  were  designed  to  charge  (in  phase- 1)  and 
discharge  (in  phase-2)  the  bridging  oxygen  vacancy  among  its  three  charge  states  +1, 
0,  -1  in  order  to  give  controlled  variation  of  Coulomb  scattering  of  the  channel 
electrons  which  limits  the  electron  mobility.  The  results  demonstrated  reversible 
degradation  in  all  the  aforelisted  transistor  and  material  parameters  as  anticipated  by 
device  and  material  physics.  In  phase- 1  of  the  SAM  experiment,  the  mobility  was 
observed  to  decrease  during  the  positive  oxide  charging  (+Qot)  by  electron  impact 
emission  from  the  E^-VeV  level  of  the  neutral  oxygen  vacancy  center,  Vq  +  e"*  -^  Vq 
4-  2e~  and  it  also  decreased  after  negative  oxide  charging  (-Qqt)  by  thermal  electron 
capture  at  77°K,  Vq  +  e"  ->  Vq.  In  phase-2,  the  mobility  (and  also  all  the  other 
parameters  listed  above)  recovered  by  low-field  SHEi  and  subsequent  thermal  capture 
at  the  previously  positive  charged  oxygen  vacancy,  Vq  -i-  e~  ^  V^.  The  mobility  also 
recovered  by  isochronal  electric  field  stimulated  electron  emission  from  the  negatively 
charged  oxygen  vacancy,  Vq  — >  Vq  +  e~.  A  novel  methodology  was  demonstrated  also 
in  this  chapter  to  separate  the  A]x^^^  and  AV^j  contributions  to  the  degradation  in  I^.^^^, 
by  using  a  two-parameter  ([i,i„  and  Wqj)  least-squares  fit  to  the  Sah-Pao  MOS  transistor 
drain-current  equation  that  included  the  bulk  or  body  charge  [55].  The  anticipated 
linear  dependence  of  the  channel  electron  mobility  on  the  strength  (density)  of 
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Coulombic  scattering  by  the  oxide  charge  is  demonstrated  experimentally  in  this 
thesis. 

The  second  part  of  this  thesis,  chapter  5,  demonstrated  a  new  highly  accurate  d.c. 
technique  for  separating  the  effects  of  areal  oxide  trap  nonuniformities  in  the  y- 
direction  (along  the  channel)  and  the  presence  of  interface  traps  on  the  distortion 
exhibited  in  the  IV  and  CV  characteristics  of  a  MOST.  The  technique  involved  a 
comparison  of  the  two  d.c.  measurements:  (1)  the  subthreshold  Iq-Vq  swing  and  (2)  the 
new  DCIV  base  recombination  current,  le-V^jg.  The  experiments  were  conducted  in 
two  phases  which  first  generated  areally  uniform  or  nonuniform  oxide  charge  Qqj  in 
the  gate  oxide  layer  and  interface  trap  Qj-p  at  the  SiOj/Si  interface.  The  oxide  charge 
was  then  neutralized  by  uniform  low-field  SHE  injection  or  EFSE  emission.  The 
initial  positive  charging  of  the  neutral  oxygen  vacancy  by  electron-impact  emission  of 
a  bound  electron  at  E^  -  7eV  energy  level,  Vg  +  e"*  — >  Vq  +  2e~,  was  shown  to  be 
areally  uniform  during  high-field  SHEi  (V^^-OV)  but  areally  nonuniform  during 
SHEi/CHEi  (Vps=4V)  and  FNTEi  under  accumulation  and  inversion  d.c.  gate  voltage 
polarities.  The  low-field  SHEi  at  77°K  during  phase- 1  was  shown  to  give  nonuniform 
negative  charging  of  the  neutral  oxygen  vacancy  by  electron  capture  into  the  shallow 
acceptor  level  at  E^  -  leV,  Vq  +  e"  -4  Vq-  The  key  assumption  made  (validated  by 
these  experiments)  in  demonstrating  this  technique  was  that  any  increase  in  the  density 
of  interface  traps  could  not  be  removed  during  the  uniform  neutralization  of  the 
trapped  oxide  charge  generated  during  phase- 1. 

Finally,  a  new  current-accelerated  CHE  stress  method  for  rapid  time-to-failure 
extrapolation  or  determination  was  demonstrated  in  chapter  6.  The  advantages  of  this 
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current-acceleration  method  over  the  traditional  voltage-accelerated  CHE  stress 
method  is  that  the  new  method  keeps  the  kinetic  energy  and  spatial  distributions  of  the 
channel  electrons  constant  to  insure  that  the  same  mechanisms  are  responsible  for  the 
change  or  degradation  of  the  transistor  device  and  materials  parameters  (listed  above). 
In  contrast,  the  tradition  voltage-acceleration  CHE  method  greatly  alters  both  the 
kinetic  energy  and  spatial  distributions  of  the  electrons  and  hence  also  changes  the 
fundamental  degradation  mechanisms  during  the  stress  tests.  Stress-time  acceleration 
factors  of  more  than  200x  are  demonstrated  for  the  new  current-accelerated  CHE  stress 
method  in  this  thesis  by  forward-biasing  the  bottom  emitter  of  the  n-BiMOST  or  the 
source  n-i-/p  junction  as  the  top-emitter  to  increase  the  flux  of  electrons  flowing  through 
the  channel.  The  saturating  behavior  of  the  density  of  the  interface  traps  generated 
near  the  drain,  exhibited  in  the  Alg-t^j^^^,  data,  gives  for  the  first  time  a  demonstration  of 
the  anticipated  asymptotic  limit  of  CHE  degradation  at  lower  voltages  and  long  stress 
times,  which  degradation  clearly  cannot  increase  indefinitely  as  has  been  tacitly 
assumed  in  the  traditional  linear  or  log-log  empirical-statistical  fit  used  by 
manufacturing  engineers  to  develop  future  generations  of  deep  sub-half-micron  silicon 
MOS  integrated  circuit  technologies. 

A  major  practical  application  of  the  new  measurement  methods  and  new  results 
demonstrated  and  presented  in  this  thesis  is  Reliability  Computer-Aided  Design 
(RCAD)  for  accurate  prediction  of  the  reliability  of  future  deep-submicron  ULSI 
technologies  under  development.  This  physics-based  new  RCAD  method  will  provide 
reliable  fundamental  degradation  rate  parameters  at  the  anticipated  low  operation 
voltages,  so  that  the  RCAD  algorithms  would  not  rely  on  the  traditionally  used 
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empirical-parameter  extraction  procedures  to  calibrate  the  CAD  tools  which  is 
expected  to  be  in  serious  error  due  to  the  saturation  or  cessation  of  the  fundamental 
degradation  mechanisms  and  pathways  at  the  low  integrated  circuit  operation  voltages 
(<3.3V  to  less  than  IV)  and  the  corresponding  low  electron  kinetic  energies.  A  major 
economics  impact  comes  from  the  more  than  lOOx  acceleration  of  the  new  current- 
acceleration  method  which  will  shorten  the  time-to-failure  tests  for  new  MOS 
technologies  being  developed,  from  the  three-to-six  months  required  in  current  factory 
practices  to  less  than  100  hours,  thereby  substantially  speeds  up  the  development  and 
production  cycle  of  new  generations  of  multi-million  transistor  silicon  MOS  integrated 
circuits. 
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